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ABSTRACT 


When  designing  an  active  sonar  homing  torpedo,  certain 
operational  torpedo  parameters  such  as  speed,  turn 
rate,  etc.  have  to  be  decided  upon.  For  a  given 
homing  torpedo,  there  must  exist  tactical  guidelines 
of  how  to  employ  the  torpedo,  i.e.  which  firing 
position  gives  the  best  chance  of  a  hit.  This  thesis 
attempts  to  gain  some  insight  into  the  detection 
process  during  the  torpedo  run,  as  well  as  getting 
some  indications  of  the  relative  importance  of  the 
different  torpedo  parameters  and  the  tactical 
situations.  A  simulation  model  was  used  in  order  to 
generate  the  data  base  for  analysis.  The  results 
stress  the  importance  of  a  good  firing  position  as 
well  as  show  how  it  is  possible  to  counter  a  bad 
firing  position  by  a  high  speed  torpedo.  They  also 
point  to  the  importance  of  having  only  cne  detection 
as  requirement  for  target  acquisition. 
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I.   INTRODUCTION 


The  following  analysis  examines  the  performance  of  a 
homing  torpedo  against  a  surface  ship.  A  homing  torpedo  is 
described  as  a  torpedo  which  is  searching/snaking  on  each 
side  of  its  main  course.  It  is  searching  for  a  target  by 
transmitting  with  its  sonar  and  listning  for  an  echo.  Ref. 
Fig.  1.  Passive  searching  torpedoes  and  homing  torpedoes 
going  in  circles  are  not  investigated  in  this  paper. 

The    torpedo's    performance   is   a   function   of   many 
variables.   These  variables  are  divided  into  two  groups; 
technical  variables;  speed,  max  torpedo  run,  sweep 
angle,  technical  detection  range,  lobe  characteristics 
and  turn  rate. 

tactical  variables  ;  firing  range,  attack  angle, 
target  speed,  type  of  target  and  tactical  detection 
range  (sonar  conditions) . 

No  attempt  is  made  to  analyze  the  first  group  of 
variables;  instead,  technical  variables  used  are  those  of 
present  technology.  We  are  assuming  a  'standard  homing 
torpedo'  based  upon  homing  torpedoes  in  operational  use 
today  [8],  This  'standard  torpedo'  assumes  conventional 
warhead  and  active  sonar  transmission  for  detection.,  and  is 
unguided. 
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The  technical  variables  (torpedo  parameters)  are  ir.  many 
ways  interrelated.  For  example,  the  maximum  detection  range 
will  determine  the  transmission  rate,  since  the  transmitted 
energy  must  have  time  to  traverse  out  to  maximum  detection 
range  and  return  as  an  echo  before  the  next  transmission,  at 
least  during  the  search  phase. 

At  the  same  time  the  torpedo  is  transmitting,  it  is 
searching  (changing  course)  for  a  target.  In  each 
transmission/  the  transmitted  energy  is  focused  within  a 
narrow  beam  (lobe)  .  During  reception,  the  echo  is  confined 
within  the  same  narrow  beam  (lobe).  Concurrently,  in  the 
time  between  two  transmissions  the  turning  rate  of  the 
torpedo  must  be  limited  to  ensure  that  the  receiving  lcbe  is 
not  outside  the  direction  from  where  an  echo  may  return. 
Thus  turn  rate  should  be  a  function  of  detection  range  and 
the  lobe  pattern. 

In  order  to  maintain  torpedo  speed,  the  number  of 
degrees  of  sweep  on  each  side  of  the  main  course  must  be 
small.  If  the  sweepangle  is  small,  however,  the  width  of 
the  possible  detection  lane  will  be  small  as  well,  and 
consequently  the  detection  probability  might  be  reduced 
during  transit.  Also,  a  high  torpedo  speed  creates  a  great 
change  in  torpedo  position  between  each  transmission.  In 
this  way  the  torpedo  may  scan  outside  a  target  in  the  sweep 
lane.  In  other  words,  the  coverage  density  cf  the  lobe  may 
be  low  as  a  result  of  the  high  movement  rate. 

As  we  recognize  the  relationship  between  torpedo 
parameters,  tactical  variables  and  torpedo  performance,  we 
know  that  frequently  within  the  naval  establishment 
decisions  have  to  be  made  with  regard  to  torpedo  parameters 
and  tactical  doctrines.  In  localizing  and  defining  these 
relationships   this    analysis    may   be   a   tool   in   this 


13 


decision pro cess. 

The  measure  of  effectiveness  by  which  different 
alternatives  will  be  judged  will  be  detection  probability, 
by  which  is  meant  the  probability  that  the  torpedo's  active 
sonar  detects  and  begins  to  track  the  target.  This 
probability  will  be  measured  by  a  digital  computer 
simulation,  construction  of  which  was  a  major  part  of  the 
author's  effort  in  writing  this  thesis. 
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II.   NATURE  OF  TH2  PROBLEM 


A.   DEFINITIONS 


Lobe  width  is  the  number  of  degrees  from  the 
centerheading  of  the  torpedo,  until  the  first  minimum  in 
transmission  intensity  is  reached.   See  Fig.  4. 

Detection  range  is  the  range  to  the  target  when 
detection  first  occurs. 

Technical  detection  range  is  the  max  detection  range 
which  is  technically  and  reasonably  possible  considering 
power  transmitted  and  lobewidth.  It  is  the  basis  for 
determining  the  transmission  rate. 

Aspect  is  the  angle  measured  from  the  positive  direction 
of  the  longitudinal  axis  of  the  target  to  a  line  joining  the 
centers  of  gravity  of  the  target  and  the  torpedo. 

Attack  angle  is  the  aspect  at  the  start  of  the  torpedo 
run. 

Sweep  angle  is  the  maximum  number  of  degrees  the  tcrpedo 
will  turn  off  the  main  course  during  search. 


All  dimensions  are  in  meter,  second,  meter  per  second, 
degree,  degree  per  second.  Speed  of  the  target  and  the 
torpedo  are,  however,  always  given  in  knots. 
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It  is  assumed  that  all  firings  are  successful,  and  the 
torpedo  will  not  deviate  from  its  ordered/calculated  course 
and  speed. 

All  firings  are  made  with  a  deflection  angle;  i.e.  the 
torpedo  is  given  a  course  to  a  predicted  hitting  point  with 
the  target. 


B.   ASSUMPTIONS 


Not  only  in  order  to  keep  the  problem  tractable,  but 
also  because  of  modern  torpedo  development,  only  surface 
targets  are  considered.  Previously  within  NATO,  torpedo 
developments  seemed  to  start  as  a  development  of  an 
anti-submarine  torpedo  with  later  modifications  in  order  to 
make  the  torpedo  dual  purpose.  However,  today  there  are 
some  indications  that  the  anti-surface  ship  requirement  is 
coming  into  the  development  early  in  the  planning  process 
[7;8;9].  The  entire  problem  is  then  kept  in  two  dimensions. 
The  vertical  axis  is  not  significant  as  we  assume 
isovelocity  condition,  and  we  assume  for  simplicity  that  we 
have  negligible  surface  effect. 

Also,  if  the  intensity  of  the  echo  is  above  detection 
threshold  level,  the  target  is  detected  with  probability 
one.   Probability  of  false  contact  is  assumed  to  be  zero. 

The  main  purpose  of  a  homing  torpedo  is  to  counter 
uncertainty  in  target  data  at  firing  and  target  maneuvering 
after  firing. 

For  simplicity  the  following  assumptions  are  made: 

the  target  remains  on  a  steady  course  after  firing, 
estimated  target  data  is  used  in  solving  the 
deflection  angle  problem 
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deflection  angle (DA)  is  given  by; 

DA  =  ARCSIN((TAM  X  SIN  (ASP)  )  /TO)  (2.1) 

TAM    =  target  estimated  speed 
ASP     =  estimated  aspect 

=  (target  estimated  course)  - 

(bearing  to  torpedo) 
TO     =  torpedo  speed. 
See  Fig..  2. 

The  difference  between  the  target  data  and  the  target 
estimated  data  are  defined  as  errors  in  the  target  data. 
These  errors  are  assumed  to  be  random  variates  and  are  given 
as; 

-   target  range  error  is  uniformly  distributed  between 

-  15  %   and  +  15  %   of  actual  target  range 

target  course  error  is  uniformly  distributed  between 

-  15  and  +  15  degrees 

target  speed  error  is  normally  distributed  with  mean 
0  and  standard  deviation  3  knots. 

These   errors   are   assumed   to  cover  errors  in  the  fire 

control   solution   at   the   time   of   firing   as  well   as 

non-radical   maneuvering   of   the   target  during  the  tcrpedo 
run. 

As  shown  in  Eg.  2.1,  estimated  range  does  not  enter  into 
the  calculation.  Estimated  range  would  only  be  used  for 
some  more  complicated  tactical  situations  as  angled  torpedo 
firing  off  the  firing  course  of  the  firing  unit.  These 
situations  are  not  covered  in  this  study. 
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Figure  2  -   TORPEDO  TRIANGLE 
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III.   PROBLEM  SOLVING  APPROACH 


At  firing,  the  initial  course  of  the  torpedo  is 
uniformly  distributed  between  minimum  course  and  maximum 
course  (main  course  +/-  a  fraction  of  sweep  angle) . 

Immediately  after  firing,  the  torpedo  starts  'snaking*. 
During  snaking,  the  torpedo  is  continously  changing  course 
left  or  right  out  to  the  given  sweepangle,  then  back  pas- 
main  course  and  out  to  sweepangle  on  the  other  side  and  so 
on.  The  torpedo  is  turning  with  the  given  turnrate.  Euring 
the  whole  process,  the  torpedo  is  also  transmitting  and 
listening.  Transmission  interval (TTIME)  is  given  by 
technical  detection  range  as; 

TTIME  =  2  x  TEDEC/1500      seconds  (3.1) 

where 

TEDEC  =s  technical  detection  range  in  meters. 
1500   =  speed  of  sound  in  salt  water,  m/sec. 

The  torpedo  run  is  conducted  in  steps.  Every  0.5 
seconds  interval,  all  positions  and  courses  are  updated. 

At  each  transmission;  the  relative  bearing  to  target, 
and  the  target  aspect  are  calculated  in  order  to  establish 
the  intensity  of  the  echo. 

When  a  detection  occurs,  the  following  data  is  stored; 
detection  range  to  the  center  of  the  target. 
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-    detection  range  to  the  nearest  part  of  the  target. 
detection  bearing  (relative)  to  the  center  of  the 
target. 

detection  bearing (relative)  to  the  nearest  part  cf 
the  target, 
target  aspect. 

In  addition  to  the  detection  probability,  the  range  at 
which  the  detection  first  occurs  is  also  of  interest. 
Therefore,  we  store  these  data  at  the  first  detection. 

However,  successive  detections  are  also  important.  As 
part  of  the  criterion  for  the  decision  of  when  to  go  from 
search-phase  to  attack-phase,  the  number  of  successive 
detections  (with  no  non-detection  between)  may  be  employed. 
In  real  life  there  is  always  a  positive  probability  of  false 
detection.  Even  if  we  are  not  addressing  the  problem  of 
false  contact  as  such,  we  can  cover  the  possibility  by 
requiring  the  torpedo  to  have  at  least  two  successive 
detections  before  going  into  attack-phase.  Accordingly,  we 
store  also  the  previously  listed  data  at  the  second 
successive  detect  ion  (two  immediately  following  detections), 
at  the  third  and  so  on,  up  to  and  including  5  successive 
detections.  This  listing  of  detections  will  give  an 
indication  of  the  decrease  in  detection  probability  if  a 
large  number  of  successive  detections  before  going  into 
attack-phase  is  required  in  order  to  decrease  the 
probability  of  false  contacts. 
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IV, 


MODEL 


A.   SEARCH 


For   simulating   the  torpedo   search,   a   Fortran   IV 
simulation  program  was  developed. 
The  program  was  divided  into; 

Main  program,  including  generation  of  statistics  and 

print  out  of  summary  after  all  the  runs  were  completed. 

Subroutine  PARMET  for  setting  tactical  situation  and 

torpedo  parameters. 

Subroutine  FIRING  which  calculates  estimated  target 

data,  and  the  deflection  angle. 

Subroutine  POSIS  which  calculates  the  torpedo  course, 

and  torpedo  and  target  positions  at  each  time  step. 

Subroutine  DETECT  which  checks  if  the  target  is 

detected  and  if  so,  store  detection  data. 
See  Fig.  3. 
See  Appendix  A  and  Appendix  B. 


21 


MAIN 
PROGRAM 


SUBROUTIN 
PARMET 


SUBROUTINE 
FIRING 


SUBROUTINE 
POSIS 


SUBROUTINE 
DETECT 


Figure  3  -   STRUCTURE  OF  COMPUTER  PROGRAM 
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B.   DETECTION  MODEL 


A  contact  occurs  when  the  acoustic  energy-pulse 
generated  at  the  transducer  and  reflected  from  the  target  as 
an  echo,  is  at  or  above  threshold  level.  In  the  following 
discussion  we  assume  that  the  contact  meets  the  tactical 
requirement,  and  accordingly  we  use  the  term  detection. 

1 •   Detection  Threshold 

Deciding  if  a  detection  occurs  is  a  function  of 
detection  threshold  (signal  to  noise  ratio),  the  range  to  the 
target,  the  target  strength  and  the  relative  bearing  to  the 
target,  given  a  level  of  radiated  intensity. 

The  detection  threshold  for  a  torpedo  is  a  function 
of  design  and  technological  sophistication  of  the  torpedo. 
Without  making  any  assumption  about  these  variables  in  the 
model,  we  start  with  a  given  technical  detection  range,  a 
•standard'  target,  and  calculate  intensity  of  echo  at  that 
range  for  target  aspect  equal  to  90  degrees  (maximum  target 
strength)  and  relative  bearing  to  the  target  equal  tc  zero 
degrees.  This  echo  intensity  is  then  the  detection 
threshold  for  every  transmission  during  a  run.  If  any  echo 
intensity  is  above  the  detection  threshold,  it  is  detected; 
if  below  the  detection  threshold,  it  is  not  detected. 

2.   Echo  Intensity 

In  calculating  echo  intensity  we  must  separately 
investigate  the   important   factors,   which   are   transducer 
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gain,      lobe      characteristic,      transmission      loss      and    target 
strength.      The   model    used    is   described    below. 


2U 


a.   Lobe  Characteristics 

The   transducer   has   a  main   lobe    and    many 

sidelobes   as   a   function   of   the  transducer's   gain   and 

relative  bearing.   Urick  [6;51-57]  discusses   some   of   the 
different   types   of  beam  pattern  (lobes) ,  and  the  following 

mathematical  model   was   developed  and   found   to   give   an 
acceptable  pattern; 


I  SIN  (x  TT)  | 
G(0)  =  G   | |  COS  (6/2)  (4.1) 


0 


x  tt   | 


where 


and 


x  =  9/9  (4.2) 

0 


G   =  maximum  gain 
0 


9   =  relative  bearing 

9   =  lobe  width. 
0 


This   model   will   produce   the   gain-pattern  as 
shown  in  fig.  4. 
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Figure  4  -   DISTRIBUTION  OF  LOBES  AND  INTENSITY 
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b.   Reduction  in  Intensity  due  to  Range 

Primarily f  the  reduction  is  due  to   two   effects; 
spherical  spreading  and  absorption. 

Spherical  spreading  is  a  known  function,  but  absorption  is 
dependent  upon  transmission  frequences,  water,  salinity  etc. 
In  order  to  simplify  the  model  and  since  spherical  spreading 
has  the  greatest  effect,  only  the  spherical  spreading  for 
reduction  in  intensity  is  considered.  This  reduction  in  a 
one  way  propagation  is  given  by; 


I  =  I  /R  (4.3) 

0 


where 


I  =  radiated  intensity  at  one  meter 
0 

I  =  intensity  at  range  R. 

R  =  range  in  meters. 


c.   Target  Strength  and  Target  Aspect 

When  the  transmitted  energy  pulse  hits  the 
target,  some  of  the  energy  is  reflected  back,  to  the 
transducer.  The  echo  intensity  is  a  function  of  the  shape 
and  dimension  of  the  target,  type  of  reflective  material  and 
aspect. 

It  should  be  noted  that  the  notion  of  target 
strength  represents  the  ratio  between  target  cross  section 
and  the  surface  of  a  sphere  of  radius  1  meter,  or  if  in  dE, 
10  times  the  log  of  this  ratio;  base  10.  In  most 
references,  the  target  strength  or  the  target  cross  section 
is  given  abeam  of  the  target,  see  [ 5 ; 97  ], [ 6  ;274  ],  without 
presenting  the  cross  section  as  a  function  of  the  aspect. 
(Jrick  [6;282,283]  gives  ,  however,  as  figures,  an  indication 


27 


of  how  the  target  strength  (in  dB)  varies  with  the  aspect. 
Cox[3;60]  states  that  it  will  vary  between  *,0  and  25  dB. 
All  measurements  in  dB  in  the  two  references  are  relative  to 
1  yard  as  unit  for  range.  Orick[ 6 ; 283-286  ]  indicates  that 
his  reference  (as  given  in  the  figures)  will  not  change  in 
any  considerable  degree  with  changes  in  frequencies  (20-60 
KHz)  or  for  different  targets  (submarines/surface  ships) . 

Assuming  a  torpedo  with  transmitting  frequency 
between  50  and  60  KHz,  we  get  a  wavelength  varying  between 
2.5  and  3.0  cm (0.025  -  0.03  meters).  As  any  reflection  from 
a  target  is  mostly  determined  by  target  form,  size,  aspect 
and  wavelength,  we  may  use  a  model  from  radar  theory  in  our 
next  step.  The  justification  for  this  use  is  that  in  radar 
theory  we  are  working  in  the  same  area  of  wavelength  and 
target  dimension  as  an  active  sonar  for  a  homing  torpedo. 

Crispin  and  Siegel  [4; 86]  give  for  target  cross 
section  a  model  for  an  ellipsoid  where  the  incident 
angle  (target  aspect)  is  a  variable.  The  relationship  is  as 
follows; 


6    = 


2   2   2 
TT-a  -b  -c 


(U-4) 


ar  b,  c  being  half  axes  of  the  ellipsoid. 
Ref.  Fig.  5. 
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Figure  5  -   MODEL  OF  TARGET  AND  TARGET  ASPECT 
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As  we  assume  that  the  transmitting  pulse  is 
always  in  the  horizontal  plane,  9  is  90  degrees,  which  gives 
us; 


_     2      2      2 
IT -a     -b    -c 


o  = 


(a    -COS    0    +    b     -SIN  $)  ' 


(4.5) 


(p    =   aspect. 

Urick      [6;275]    gives   for   target   section    a    model 

for   abeam  or    ahead   cases,    which    is; 

2 
t    =    d/  (4TT)    =     (b  -c/2-a)     , 

identical   with   Eg.   4.5.    Note   that   Eg.    4.5    is   an 

expression  for  the  target  cross  section. 

Haslett  [5; 139]  gives  for  the  target  cross  section   a   model 

for   both   ahead  and  abeam  cases.   His  model  equals  Eg.   4.5 

2 
times   a   factor   R  ,   where   R   is   acoustic   reflectivity 

coefficient  (per  cent)  =94. 

The  advantage  of  using  Eg.  4.5  is  that  it  gives  the  target 
cross  area  as  a  continous  function  of  the  target  aspect. 
For  our  model  we  will  only  use  the  lower  part  of  the 
ellipsoid  to  simulate  the  ship  hull  below  the  water  line. 

Combining  Eg.  4.5  and  acoustic  reflectivity 
coefficient  we  get  the  following  model  for  the  target  cross 
section; 


0  = 


2   2   2   2 
Tf-a  -b  -c  -R 


(a  -COS  p  +  b  -SIN  p)  4 


(4.6) 


With    reference    to    Urick's    figures   [6;283]   where 
the    pattern    of    the    target    strength    is    given    as   a    function   of 
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aspect  in  Fig.  9.13,  and  reproduced  in  this  analysis  as 
Fig.  6. a,  we  still  have  not  obtained  a  model  which  gives 
the  same  type  of  pattern.  By  applying  the  following  scaling 
factor  to  Eg.   4.6   we  have  approximated  his  information: 


2 
D  =  (0.251635-1)   -  0.18555-(p  +  0 .036  5  •  SIN  (3  •  (<p  ♦  0.17453)) 


2  -1 

+  0.015-p  -SIN(9- (p/2)  )  (4.7) 


We  then  have  as  the  target  cross  section  in  oar  model  the 
following   expression; 

0   =   tf    x   U  (4.8) 

where  0  and  U  are  as  previously  shown. 

Fig.  6. a  and  Fig  6.b  shows  an  'ideal1  pattern  and  a  model 
pattern.  The  figures  given  by  (Jrick  are  for  1  yard  as 
reference  distance,  but  have  been  converted  to  1  meter 
reference  distance  in  Fig.  6.  To  go  from  d3(yard)  to 
dB  (meter) ,  we  subtract  0.78  dB.  The  dB  as  given  in  this 
analysis  are  all  with  1  meter  as  reference  distance. 
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Figure  6  -   TARGET  STRENGTH 
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The    active    sonar    equation    is; 


P   = 


p    -G    -O  -G     •  X 
0      t   "      r 

3 — IT — 

(4  If)      -R 


Watts 


(4.9) 


P   =  power  received 

P   =  power  transmitted 
0 

G   =  gain  transmitting,  ref  Eg.  (4.1) 

d  =  target  cross  section,  ref  Eg.  (4.8) 

G   =  gain  receiving,  ref  Eq.  (4. 1) 

r 

A  =  wavelength  in  meters. 

R   =  range  to  target  in  meters. 

This   may   be  rewritten  into  an  expression  of  power  received 
as  a  function  of  the  variables  of  the  different  terms; 


SIN  (X  -If)  SIN  (X  -If) 

t        2   2   2   2        r 
|  -a  -b  -c  -R   | 

■~x™rir        |  -    i  — *~ir 

t  r 


P  =  K 


R  .(a  -COS  p  +  b  -SIN  <p) 


(U.10) 


K  =  the  product  of  all  the  constants  in  the  tercrs. 

For   more   detailed   discussion  about    gain, 

transmission  loss  and  reflection  (target  cross  section) ,  see 
C 1  ;  110-1 11]  and  [ 6  ;29, 94, 263  ]. 

We   can   now  calculate  the  minimum  power  level  for  detection 
by  setting: 
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R   =  technical  detection  range 
X   =0  degree 

X   =0  degree 

r 

<p      =90    degrees 
and   we   get 


2      2      2      2 
a    -b     -c    -R     -0 
P  =    K (4.11) 

min 

4         2    2 
B     -(b   ) 

t 


and  by  substituting  for  0 


2   2   2   2 
a  -b  -c  -R  -3.08657 

P     =  K (4.12) 

min 

4    2  2 
R  .(b  ) 
t 

a,   b   and   c   are   the   dimension  of  the  target  used  in  the 

model. 

He   assume   a   'standard*  target,  length  100  meters,  beam  15 

meters  and  draught  4  meters,  i.e. 

a  =  100 

b  =  15 

c  =  4. 


There  will  be  a  detection  if  p/p    >   1.    Note 

min 

that   P/P    does  not  depend  on  K,  into  which  radiated  power 
min 

and  transducer  gains   have   been   included.    The   technical 

detection   range  R   is  an  operationally  meaningful  surrogate 

t 
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for  these  parameters.   The  ratio  P/P     will  be   called   the 

min 

"intensity  fraction". 
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3 .   Detection  Rule 

Any  intensity-fraction  calculated  during  a 
transmission  which  is  greater  than  1  is  a  detection. 
However,  to  improve  the  model  at  close  ranges,  the  following 
modification  has  been  made  for  gain  variation  due  to 
relative  bearing. 

At  close  ranges,  the  relative  bearing  to  target  can 
alter  considerably  from  bow  to  stern.  Therefore,  the 
intensity  in  the  pulse  will  differ  along  the  target.  To 
average  this  intensity  both  for  the  radiated  pulse  and  for 
the  echo,  the  model  calculates  relative  bearing  to  the 
target  bow,  center  and  stern,  calculates  the  corresponding 
gain  factor  for  each  bearing,  and  finds  the  arithmetic  mean 
of  these  gain  factors.  These  two  average 
gainfactors  (transmitting  and  receiving)  are  then  used  in  the 
calculation  of  echo  intensity. 

The  model  does  not  recognize  a  detection  unless  the 
tactical  situation  makes  it  possible  to  maintain  contact 
with  the  target  for  some  time.  To  be  precise,  the  following 
conditions  must  be  present; 

-  torpedo  turn  rate  higher  than  bearing  rate 

-  closing  speed  must  be  positive. 

-  target  must  have  2  knots  doppler. 


36 


7,   PRESENTATION  OF  DATA 


A.   STOCHASTIC  ELEMENTS 


In   the  previous  description  of  the  model,  the  following 
input  values  are  stochastic; 
error  in  target  speed 
error  in  target  course 
error  in  target  range 
-   initial  torpedo  course (not  main  course). 


The  primary  stochastic  effects  on  the  torpedo 
performance  are  identified  as  errors  in  target  speed  and 
course,  since  these  two  variables  are  the  only  stochastic 
ones  used  in  computing  the  torpedo  main  course.  The  first 
problem  to  be  solved  was  then  how  to  design  the  run  series 
in  order  to  reduce  variance  in  result  at  the  same  time  as 
keeping  the  result  unbiased. 

It  was  found  that  instead  of  using  a  complete  randomized 
design  (random  variates) ;  we  could  deterministically  section 
the  probability  range  0.0  -  1.0  for  the  two  important  random 
variables,  using  the  inverse  probability  transformation  to 
get  variates,  and  then  run  the  number  of  runs  required  to 
cover  all  combinations  of  variates. 


Some  preliminary  simulation   runs   were   done   in   three 
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versions;  complete  randomized  and  independent;  with 
antithetic  reduction  technique  (sectioning) ;  and  the 
previously  described  procedure.  The  number  of  runs  needed 
in  order  to  keep  the  variance  low  for  the  result  was 
considerably  higher  for  the  first  two  versions. 
Accordingly,  we  selected  the  previously  described  procedure. 
It  was  found  that  a  series  of  150  runs  was  sufficient  in 
order  to  give  a  reasonable  accuracy  in  detection  probability 
and  at  the  same  time  keeping  the  total  CPU  time  for  a  series 
of  runs  acceptably  low.  The  150  run  series  was  established 
by  dividing  the  range  of  probability  of  target  speed  errors 
into  15  equally  spaced  sections;  and  the  rang€  of 
probability  of  target  course  errors  into  10  equally  spaced 
sections.  Each  section  boundary  point  was  by  inverse 
probability  transformation  converted  into  a  variate. 
Bearing  in  mind  that  speed  errors  are  normally  distributed 
and  course  errors  are  uniformly  distributed,  all 
combinations  of  target  data  error  are  plotted  in  Fig.  7. 
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Figure  7  -   DISTRIBUTION  OF  ERROR  IN  TARGET  DATA 
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B.   TYPE  OF  PRINTOUT  OF  DATA  AND  RESULT 


Each  series  of  150  runs  produces  a  printout  as  shown. 
The  heading  cf  the  printout  gives  the  tactical  situation  and 
the  torpedo  parameters  in  the  given  run  series.  Also  ,  the 
printout  gives  the  sweep  lane,  which  is  the  width  of  the 
lane  where  the  torpedo  has  swept  through  by  its  sonar  lobe. 
The  coverage  ratio  gives  an  indication  of  the  fraction  of 
the  lobe,  which  is  covered  twice;  i.e.  how  much  the  lobe  is 
being  offset  from  its  previous  position  by  change  in  the 
torpedo  course.  The  question  of  offsetting  the  sonar  lobe, 
about  which  information  is  given  in  the  printout,  is 
discussed  later  in  Ch.  VI. 
Ref.  Fig.  8. 

For  each  run,  the  following  are  output:  target  data, 
torpedo  deflection  angle,  torpedo  main  course,  target  and 
torpedo  grid  position  at  end  of  run,  duration  of  torpedo  run 
and  length  of  torpedo  run. 

After   all   runs   in   a  series   are  completed,  a  summary  is 

given. 

Ref.  Fig.  9. 

The   summary   gives   detection   probability   for   a    single 

detection,   2   successive   detections,   up   to   5  successive 

detections.   Also  mean  detection  range,   standard   deviation 

of   detection   range,   mean   aspect,   mean  detection  bearing 

relative  to  center  bearing  of  sonar   lobe   and  relative   to 

main  course  are  given. 


Lastly,  the  detection  range,  the  relative  bearing  to  the 
center  of  the  target  and  to  the  closest  part  of  the  target 
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Figure    9    -       EXAMPLE    OF    PRINTOUT    SUMMARY 
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at  detection  (relative  to  present  torpedo  course)  ,  and  the 
target  aspect  at  detection  are  printed  for  each  run  for  a 
single  detection,  2  successive  and  3  successive   detections. 

It   also   should   be  noted  that  it  is  possible  to  get  a 
more  detailed  printout  for  each  run  by  setting  IPRINT  =  0  in 
the  main  program  (main  program  statement  035) . 
Ref.  Appendix  D.   for  example  of  detailed  run  printout. 


From    the   printout  data,   it   is   possible   to   study 

different  aspects  of  the  detection  process   as   well   as   to 

generate   distributions  of  detection  range,  aspect,  bearing 
etc. 
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VI.   PARAMETRIC  TORPEDO  ANALYSIS 


A.   OBJECTIVES 


The  following  approach  was  used: 
The  torpedo  speed,  the  technical   detection   range   and   the 
lobe  width  were  assunred  to  characterize  a  torpedo  type. 
Within  the  type,  it  was  possible  to  change  the  turn  rate  and 
the  sweep  angle. 

A  tactical  situation  was   characterized   by   the   attack 
angle,  the  target  speed  and  the  firing  range. 


The  following  questions  were  investigated: 

Can  a  torpedo  be  improved  by  offsetting  its  sonar  lobe 

from  the  torpedo  heading  ? 

Rephrased;  it  may  be  asked,  is  the  sonar  lobe  searching 

in  the  right  direction  (most  likely  area)  by 

pointing  straight  ahead  along  the  torpedo  course  ? 

How  do  turn  rate  and  sweep  angle  affect  a  torpedo's 

MOE  ? 

How  are  the  different  torpedo  types  related  to  each 

other  with  regard  to  detection  probability  (MOE)  ? 

In  the  analysis,  we  started  with  a   reasonable   tactical 
situation;    target   speed   18   knots,   range   3000   meters, 
technical  detection  range  750  meters.   Initially,  we  changed 
the  attack  angles. 
With  regard  to  torpedoes,  we  started   with   three   types   of 
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torpedoes;  24  knots,  32  knots  and  40  knots;  all  with  20 
degree  lobe  width,  6  degree  per  second  turn  rate  and  30 
degree  sweep  angle. 


B.   OFFSETTING  SONAR  LOBE 


The  hypothesis  was  that  when  a  torpedo   is   fired   on   a 
deflection   angle   course,   the   sonar   lobe   should  be  most 
effective  if  it  scans  across  the  bearing  to  the  target. 
Or,   the   sonar   lobe   should  be  offset  equal  to  deflection 
angle  (DA).   Ref.  Fig.  10. 


It   was  found  that  offsetting  had  a  positive  effect  when 
attacking  from  ahead  of  target. 
Ref.  Fig.  11. a.  and  b. 

But  from  about  30  degree  to  about  110  degree  attack  angle 
the  effect  was  negative.  If  more  than  110  degree  attack 
angle,  there  was  no  effect. 

In  analyzing  the  fraction  of  offsetting,  we  analyzed  the 
case  of  30  degree  and  60  degree  attack  angle.  There  seemed 
to  be  no  effect  from  0.0  to  0.5  x  DA;  if  more  than  0.5  x  DA 
there  was  a  decreasing  efficiency. 

This  was  found  for  2  types  of  torpedoes  (32  and  40  knots;  20 
degree  lobe  width)  at  2  different  sets  of  turn  rates  and 
sweep  angles. 


This  conclusion  applies  for  both  single  detection  and 
multiple  successive  detections;  however,  the  magnitude  of 
the  effect  is  changing  as  we  look  on  different  number  of 
successive   detections.    The   conclusion   was  that  there  is 
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little  to  be  gained  by  offsetting  the  sonar  lobe,  and  the 
sonar  lobe  was  therefore  not  offset  in  subsequent 
investigations. 


U6 


Offset 
sonar    lobe 


Torpedo 


DA  -    Deflection     angle 


Figure  10  -   OFFSET  SONAR  LOBE 
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Figura  11  -   EFFECT  OF  OFFSETTING  SONAR  LOBE 
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Figure    11. b.    -    EFFECT    OF    OFFSETTING    SONAR    LOBE 
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C.   EFFECT  OF  TURN  RATE 


The  effect  of  turn  rate  was  investigated  in  the  range  3 
to  21  degrees  per  second  in  steps  of  3.  For  both  types  of 
torpedoes  the  model  showed  an  increase  in  MOE  as  turn  rate 
was  increased.  The  MOE  leveled  off  as  turn  rate  was 
approaching  15  -  20  degrees  per  second. 

The  reason  may  be  due  to  the  1  second  transmission 
interval  and  the  20  degree  lobe  width,  which  indicates  that 
the  torpedo  should  be  turned  at  a  turn  rate  equal  to  lobe 
width  divided  by  transmission  interval  for  maximum  MOE. 
However,  as  the  number  of  successive  detections  required  is 
increased,  we  get  maximum  MOE  at  lower  turn  rates. 
Fig.  12  shows  the  change  in  MOE  with  turn  rate  for  30  and  60 
degrees  attack  angles. 

From  Fig.  15  where  different  combinations  of  turn  rates  and 
sweep  angles  are  plotted  versus  MOE,  we  see  that  the  effect 
is  negligible  from  about  60-80  degrees  to  180  degrees  attack 
angle. 


A  6  degrees  per  second  turn  rate  is  compared  with  what 
may  be  termed  an  'optimal'  turn  rate  in  Fig.  13.  The 
'optimal'  turn  rates  were  established  by  the  general  trend 
from  Fig.  12  and  Table  II . a  and  II. b. 

The  following  turn  rates  were  identified  as  'optimal'; 
15  degrees  per  second  for  the  32  knots  torpedo 
18  degrees  per  second  for  the  40  knots  torpedo. 
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Tactical     Situation 
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Range 

TA     Speed 

Det.  range 


3000   m 
18    Knots 
750    m 


Torpedo 

parameters: 

TO    Speed 

32  Knots 

Sweep  angle 
30° 

Lobe    widtn 
20° 


o      One   detection 

a       Two   successive    detections 

A      Three   successive   detections 


21    Turn    rate    °/s 


30°  Attack   angle 
60° 


TO   Speed 

40  Knots 

Sweep   angle 
30° 

Lobe    widtn 
2  0° 


12 


15  18         21     Turn    rate   °'s 


Figure  12  -   EFFECT  OF  TURN  RATE 
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Tactical    Situation: 


Range 
TA    Speed 

Det.  nange 


3000   m 
18   Knots 

750    m 


n 
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Q_ 


c 
o 

U 
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a 


Torpedo 
parameters 

TO    Speed 

32  Knots 

Sweep  angle 
30° 

Lobe    width 
20° 

Turn    rate 


57s 

15°/s 


0°         30°        60°        90°        120°       150°      180°  Attack   angle 
°     One   detection  A     Three  successive   detections 

a     Two  successive   detections 


Torpedo 
parameters; 

TO    Speed 

40  Knots 

Sweeo  angle 
30° 

Lobe    width 
20° 


Turn    rate 
6°/s 


0 


18°/s 


30         60u       90°        120°       150°      180u  Attack  angle 


Figure  13  -   COMPARISON  OF  TORPEDOES  WITH  DIFFERENT  TORN 

RATES 


55 


We  see  here  in  Fig.  13  a  considerable  increase  in  MOE 
with  increase  in  turn  rate  for  attack  angles  less  than  60  - 
80  degrees  for  single  detection;  an  consistent  improvement 
for  2  successive  detections  in  the  same  area;  but  no  change 
or  a  slight  detoriation  for  3  successive  detections. 
It  is  quite  obvious  that  a  torpedo  which  requires  only  a 
single  detection  as  requirement  for  attack  has  a 
considerably  better  MOE,  and  a  considerably  higher  potential 
for  improvement  by  changes  in  turn  rate,  than  a  torpedo 
which  requires  more  successive  detections  for  classifying  a 
contact  as  a  target. 
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Table  II. b.  -  VARIATION  IN  TORPEDO  TURN  RATE 
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D.   EFFECT  OF  SWEEP  ANGLE 


From  preliminary  simulation  runs,  it  was  found  that  from 
90  degrees  (inclusive)  to  180  degrees  attack  angle  the 
effect  of  the  sweep  angle  was  neglegible.  The  analysis  was 
therefore  done  from  20  to  50  degrees  sweep  angle  only  for  30 
and  60  degrees  attack:  angle  for  both  the  32  and  the  40  knots 
torpedo. 
The  result  is  shown  in  Fig.  14. 


For  the  32  knot  torpedo  we  get  an  increase  from  30  to  40 
degrees  for  both  attack  angles.  From  40  to  50  degrees,  MOS 
either  levels  off  or  decrease  slowly.  As  a  conclusion,  we 
established  40  degrees  sweep  angle  as  the  'optimal'  value. 
For  the  40  knot  torpedo,  the  MOE  was  fairly  steady  over  the 
whole  range  for  30  degrees  attack  angle.  For  60  degrees 
attack  angle,  there  was  a  peak  at  30  degrees  sweep  angle, 
which  indicated  that  30  degrees  was  the  optimal  value. 


The  reason  for  the  different  sweep  angles  for  the  two 
torpedo  types  (Note;  both  have  6  degrees  per  second  turn 
rate)  may  be  due  to  the  time  it  takes  to  reach  the  target. 
The  shorter  time, the  less  area  on  each  side  of  the  main 
course  is  needed  to  be  covered  in  order  to  detect  a  target; 
i.e.  a  40  knot  torpedo  needs  only  a  30  dergree  sweep  angle, 
a  32  knot  torpedo  needs  40  degree  sweep  angle. 
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Tac  ti  cal     Situation  . 


Torpedo     Parameters 


Range 

3000  m 

TA    Speed 

18    Knots 

Det.  range 

750   m 

Lobe   width 

20° 

Turn   rate 

6°/s 

TO    Speed 
32  Knots 


TO    Speed 
40  Knots 


3  0°        40°        50° 
Sweep   angle 


o 

a 


One   detection 
Two   successive    detections 
Three   successive    detections 
30°    Attack   angle 


60s 


Figure  14  -   EFFECT  OF  SWEEP  ANGLE 


60 


The  reason  why  we  get  a  peak  and  then  a  reduction  in  HOE 
as  we  increase  sweep  angle  is  supposedly  due  to  a  sharp 
decrease  in  speed  along  the  main  course  as  sweep  angle  is 
approaching  60  degrees. 

As  example,  for  a  40  knots  torpedo  the  model  gave  35  knots 
along  main  course  for  50  degrees  sweep  angle  as  compared 
with  33.6  knots  for  20  degrees  sweep  angle.  For  a  slower 
torpedo,  the  effect  on  MOE  may  be  considerable  due  to  less 
speed  advantage  relative  to  the  target. 
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Table    III    -       VARIATION    IN    SWEEP    ANGLE 


E.   EFFECT  OF  BOTH  SWEEP  ANGLE  AND  TORN  RATE 


In  the  previous  discussion  we  changed  either  sweep  angle 
or  turn  rate  for  both  types  of  torpedo  while  we  kept  the 
other  variables  constant. 


In  plotting  MOE  for  initial  torpedo  (32  knots)  value, 
optimal  value  for  sweep  angle,  optimal  value  for  turn  rate 
and  the  'optimal'  torpedo  (having  both  the  'optimal'  turn 
rate  and  sweep  angle),  we  get  Fig.  15. 

Observe  how  the  MOE  changes  as  we  apply  the  individual 
•optimal'  values,  and  the  MOE  obtained  by  applying  both  the 
'optimal'  values. 

At  this  stage,  no  trials  were  made  in  crder  to  further 
increase  MOE  by  changing  sweep  angle  or  turn  rate  from  these 
values. 

One  essential  feature  is  that  virtually  none  of  the 
variables  so  far  have  had  any  effect  on  MOE  for  larger 
attack  angles  than  60-90  degrees. 


The  relative  difference  in  MOE  between  the  32  and  the  40 
knots  torpedo  types  is  shown  in  Fig.  16.  Both  torpedoes  are 
optimal  in  the  sense  that  the  best  values  for  turn  rate  and 
sweep  angle  have  been  chooser,  for  that  specific  speed. 
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Tactical    Situation 

Range  3000  m 

TA    Speed    18  Knots 
Det.  range    750  m 


c 
o 


Torpedo    parameters 

TO    Speed  32   Knots 

Lobe  width         20° 


One 
detection 
only 


30°        60°       90°       120°       150°      180°  Attack  angle 


Two 

successive 
detections 


0 


30°        60°        90°       120°       150°      180°  Attack  angl  e 


Figure  15  -   COMPARISON  OF  DIFFERENT  MODIFICATION  OF  A 

TORPEDO 


64 


Tactical    Situation 


Range  3000    m 

TA     Speed     13    Knots 
Det.  range     750    m 


30°        60°        90°        120°       150°       180" 

Attack    angle 
o    One   detection    only 
°    Two    successive     detections 
A    j hree    » 


TO    Speed     40  Knots 
Turn    rate      1  8  7s 
Sweep   angle    30° 


32  Knots 
15°/s 

40° 


Figure  16  -   COMPARISON  OF  TWO  DIFFERENT  TORPEDOES 
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It  is  obvious  that  the  main  differences  are  for  large 
attack  angles;  more  than  60-80  degrees. 

Especially  if  the  acquisition  requirement  is  one  detection 
only,  however,  a  32  knots  torpedo  is  slightly  better  up  to 
60  degrees  attack  angle.  This  improved  MOE  for  the  slower 
torpedo  may  be  explained  by  a  better  balance  between  the 
time  to  the  target  and  the  total  relative  speed.  A  toe  high 
relative  speed  may  prohibit  the  torpedo  from  getting  the 
target  within  its  sonar  lobe  before  the  target  is  passed. 

Generally,  however,  the  higher  speed  torpedo  is 
superior,  especially  for  larger  attack  angles  (120  degrees 
and  more) ;  this  can  be  explained  by  the  shorter  time  to  the 
target. 


EFFECT  OP  LOBE  WIDTH 


The  effect  of  changing  lobe  width  while  maintaining 
detection  range  is  shown  in  Fig.  17.  It  should  be  noted 
that  we  initially  started  the  simulation  with  an  'optimal' 
torpedo  with  20  degrees  lobe  width.  When  we  ran  the 
simulation  series  for  10  degrees  and  30  degrees  lobe  width, 
we  did  not  change  the  other  torpedo  parameters  in  order  to 
make  the  torpedo  'optimal'  for  the  new  lobe  width.  If  we 
had  carried  through  this  optimization  process,  we  might  have 
expected  an  increase  in  the  result  for  10  and  30  degrees 
lobe  width.  The  torpedo  parameter  in  question  would  most 
likely  be  turn  rate,  ref  discussion  previously  on  page  61. 

The  interesting  points  from  Fig.  17  are; 

a         10        degrees        lobe        width  torpedo  with  a 

one-detection-only      acquisition      requirement   is   as   good 


66 


as  a  20  degrees  lobe  width  torpedo  with  a 
two-successice-detection  requirement.  This  should 
indicate  what  we  have  to  pay  in  additional  power 
transmitted  when  acquisition  requirement  is  high.  Or, 
where  to  invest  research  resources;  in  transducer  or  in 
echo  filtering. 

the  equally  shaped  curves  for  increasing  lobe  width. 
However,  we  also  observe  an  increasing  difference  in 
MOE  between  the  curves  as  attack  angle  is  decreasing, 
the  importance  of  the  correct  balance  between  turn 
rate  and  lobe  width  for  successive  detections.  We 
observe  for  a  small  aspect  target  how  MOE  decreases 
drastically  when  we  reduce  lobe  width  from  20  degrees 
to  10  degrees  and  maintain  turn  rate  and  require  two 
successive  detections. 
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Tactical     Situ  a;  ion 

"A    Speed  18  Knots 

Dei   range      ^50  m 
Range  3000  m 


Torpedo    Parameter 

TO    Speed  4  0  Knots 

Sweep   ancle       30° 


Tunn    ra:e 


13  7: 


A     10°    Lobe    width 

a     20° 

o     30° 


One  detection   only 

Two  successive  detections 


Figure  17  -   EFFECT  OF  LOBE  WIDTH 
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EFFECT  OF  DETECTION  RANGE 


The  detection  range  is  a  function  of  the  design  of  the 
active  sonar  in  the  torpedo  as  well  as  sonar  condition  at 
the  time  of  the  torpedo  firing.  The  detection  range  as  a 
function  of  the  design  of  the  active  sonar  is  termed 
technical  detection  range.  The  detection  range  as  a 
function  of  both  the  design  and  the  sonar  conditions  is 
termed  tactical  detection  range,  or  just  detection  range. 
In  analyzing  the  detection  probability  as  a  function  of 
detection  range,  we  assumed  optimal  sonar  conditions  by 
equal  technical  detection  range  with  detection  range. 

Detection  range  was  varied  in  discrete  steps:  375  -  750 
-  1125  -  150C  meters. 

Figs.  18. a.  and  b.  indicate  that  detection  probability 
is  a  linear  function  of  the  detection  range  up  to  a 
detection  probability  of  0.8  -0.9  for  one  detection.  From 
the  model,  it  may  be  justifiable  to  approximate  the 
detection  probability  as  a  linear  function  from  375  m  to 
1125  m  detection  range. 

From  the  model  and  the  given  assumptions,  there  is 
little  usefulness  in  a  homing  torpedo  with  less  than  300  m 
detection  range. 

The  same  situation  is  shown  in  Figs.  19. a.  and  b.  in 
another  cut  of  the  response  surface.  We  see  here  how 
consistantly  the  MOE  has  decreased  over  the  whole  range  of 
attack  angles  when  going  from  1500  m  to  375  m  detection 
range. 
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Tactical     Situation 

TA    Speed     18  Knots 
Range  3000  m 


o       0    Attack   angle 
a    60° 

o 


Torpedo    Parameters 

TO    Speed  32  Knots 

Lobe  width  20° 
Sweep  angle  40° 
Turn    rate  15°/s 


r-      A  120 


375 


750 


1125  1500    m 

Det  range 


Figure    18    -      EFFECT    OF    DETECTION    RANGE 
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Tactical     Situation; 
TA    Speed     18  Knots 
Range  3000  m 


o         0°  Attack    angle 

a      60 

A    120 


Torpedo    Parameters. 
TO    Speed         40  Knots 
Lobe  widtn      20  ° 
Sweep   angle    30° 
Turn   rate  18  °/s 


One   detection 
Two   successive 
detections 


375 


750 


1125  1500  m 

Det.  range 


Figure    18. b.    -    EFFECT    OF    DETECTION    RANGE 
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Otherwise  the  picture  in  Figs.  19. a.  and  b.  is  as  in 
previous  similar  figures;  a  marked  decrease  in  MOE  with 
attack  angles  more  than  60  -  90  degrees,  and  with  the  faster 
torpedo  superior  over  most  of  the  range.  It  should, 
however,  be  noted  that  for  longer  detection  ranges  we  get 
maximum  MOE  at  0  degree  attack  angle  for  both  torpedo  types. 
This  effect  is  reduced  when  we  require  two  successive 
detections  for  acquisition. 

Se  also  experienced  a  considerable  decrease  in  MOE  for 
longer  detection  ranges  when  requiring  two  successive 
detections  instead  of  one.  It  seems  obvious  that  this 
reduction  is  due  to  a  larger  lateral  movement  at  the  extreme 
range.  As  noted  previously,  we  increase  the  transmission 
interval  (increase  interval  in  order  to  allcw  time  for  echo 
to  return)  when  the  detection  range  is  increased.  Keeping 
the  same  turn  rate,  the  sonar  lobe  will  turn  a  larger  angle 
between  each  transmission,  which  can  have  a  deteriorating 
effect  on  MOE  for  more  than  a  one-detection-cnly  acquisition 
requirement. 
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Torpedo    Parameters. 

TO    Speed        32     40  Knots 
Lobe  width     20    20  ° 
Sweep    angle40    30° 
Turn    rate         15     18  a/s 

I         i =►  - 


Tactical    Situation. 
TA    Speed      18  Knots 
Range  3000   m 


One   detection 


o  1500  m  Det.  range 

6    T125  m 

a     750  m 

A    37.5  m 


0°  30°        60°        90°        120°      150°      180° 

Attack    angle 


Figure  19  - 


COMPARISION  OF  TWO  TORPEDOES  WITH  CHANGE  IN 
DETECTION  RANGE 
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Torpedo    Parameters: 

TO    Speed         32    40  Knots 
Lobe  width      20    20  ° 
Sweep   angle  40    30  ° 
Turn    rate  15     1  8  °/s 


Tactical     Situation 
TA    Speed      18  Knots 
Range  3000  m 


Two   successive 
detections 


o  1500  m  Det.  range 

o  1 125  m 

a  750  m 

A  375  m 


0°  30°        60°        90°        120°       150°      180° 

Attack    angle 


Figure  19. b.  -  COMPARISION  OF  TWO  TOBPEEOES 
WITH  CHANGE  IN  DETECTION  RANGE 
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H.   COMBINED  EFFECT  OF  LOBE  WIDTH  AND  DETECTION  RANGE 


The-  following   approximate   relationships  exist  between 
lobe  width,  detection  range  and  sonar  power: 


,2 
P   G  -6    G   -A 
0  t  ~   r 
P  = Watts  (4.9) 

(4  1T)   -R 


where 

G   =  G   =  (4'1t  /W)  (6.  1) 

t    r 

2 

w  =  L 

L    =    2    x   lobe    width. 
Ref   [  1;49]. 

w  is  defined  as  solid  angle.  The  given  equation  is  valid 
for  small  lobe  width  only.  For  larger  lobe  width  the  exact 
relationship    is; 

w    =    2-1f-(1    -    COS    1)  (6.2) 

1  =  lobe  width. 
The  approximate   relationship   is  close   enough   up   to   60 
degrees  lobe  width. 

By  substituting  the  approximate   relationship   intc   Eg. 

4 
4.9   ,   we   get   a   reduction  of  L   in  receiving  echo  due  to 

change  in  lobe  width, 
or 
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4 
(L*R)   =  constant,  C3.3) 

which   combine   range   and   lobe   width,   and   implies   that 

detection  range  is  inverse  proportional  to   lobe   width   for 

constant  power  transmitted. 

It  is  therefore  possible  to  plot  this  function  for   constant 

power   transmitted,   and  use  this  as  a  prediction  of  hew  MOE 

may  change  with  change  |n  these  two  torpedo  parameters  <lobe 

width  and  detection  range) . 

This  is  done  in  Fig.  20;  and  indicated  by   the   dashed   line 

going   through   20   degrees   lobe   width   and  750  m  detetion 

range. 

We  then  ran  some  simulation  series  in  order  to  generate 
data  points  from  the  model.  The  data  points  gave  the  MOE  , 
and  by  fitting  curves  we  were  able  to  get  some  indication  of 
the  relationship  between  the  lobe  width  and  the  detection 
range  as  given  by  the  model. 

The  application  could  be  as  follows; 
For  a  given  torpedo  with  lobe  width  20  degrees  and  a 
detection  range  of  750  m,  we  ask  the  question,  can  MOE  be 
increased  without  increasing  power  transmitted  ? 
The  dashed  curve  through  the  point (20  degrees,  750  m)  is  a 
constant  power  curve,  and  by  following  the  curve  we  observe 
how  MOE  is  changing. 

From  the  figure,  it  is  obvious  that  a  narrower  lobe  and  a 
longer  detection  range  gives  a  better  result.  But  we  also 
observe  the  assymptotical  feature  of  the  curves.  We  reach  a 
point  where  the  constant  power  curve  and  the  constant  MOE 
curve  are  parallel. 

However,  it  should  be  born  in  mind  that  the  theoretical 
relationship  between  lobe  width  and  detection  range  is  an 
approximation  which  does  not  account  for  absorption-effect 
or   surface-effect.    This   implies   that  the  constant  power 
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curve  in  real  life  will  be  lowered.   Only   a 
analysis  can  say  how  much. 


more   detailed 
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Tactical    Situation 

TA    Speed  18  Knots 

Range  3000  m 


One   detection    only 


Torpedo     Parameter 

TO    Speed  40  Knots 

Sweep   angle  30  ° 

Turn   rate  1  8  °/s 


Constant    detection    probability 
Constant    power    transmitted 


20°  30° 

Lobe    widtn 


Figure    20    - 


VARIATION    IN    EFFECTIVENESS     AS    A    FUNCTION    OF 
LOBE    WIDTH    AND    DETECTION    RANGE 
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I.   EFFECT  OF  FIRING  RANGE 


The  most  important  factor  in  achieving  high  detection 
probability  is  the  difference  between  estimated  target 
position  and  actual  target  position  at  the  time  when  the 
torpedo  is  in  position  to  detect.  The  effect  on  the 
detection  probability  is  mainly  due  to  the  time  the  tcrpedo 
takes  to  reach  within  detection  range  of  target  and  the 
speed/course  errors  in  target  data. 

As  we  increased  the   firing   range,   we   experienced   as 
anticipated   a   degradation   in   MOE.    This  degradation  was 
experienced  for  both  the  32  and  the  40  knots  torpedo. 
The  variation  in  firing  ranges  were  at  the  following  values: 
1500  -  3000  -  5000  -  7000  meters. 

Fig.  21. a.  and  b.  shows  consistently  the  importance  of 
short  firing  ranges.  This  applies  to  both  one  detection  and 
two  successive  detections. 

Fig.  22. a.  and  b.  shows  an  additional  advantage  with  short 
firing  ranges;  a  considerable  improvement  at  firing  with 
small  aspect  (attack  angle),  less  than  30  degrees.  Again 
this  applies  for  both  types  of  torpedoes. 

Also,  we  get  an  indication  that  at  short  ranges,  about  1500 
meters,  there  is  no  significant  difference  in  MOE  of  the  two 
types  of  torpedoes  up  to  an  attack  angle  of  90  degrees. 
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Tactical 
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.Q 

O 

L_ 


c 
o 

-t-» 
u 

Q 


Situation: 
TA  Speed 
Det.  range 


18  Knots 
750   m 


Torpedo    Parameters: 
TO    Speed         32  Knots 
Lobe    width       20° 
Sweep   angle    40° 
Turn  rate~         15°/s 


1500  3000  5000  7000  m 

Firing    range 

Figure  21  -   EFFECT  OF  FIFING  RANGE 
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Tactical    Situation  Torpedo     Parameters: 

TA  Speed    18  Knots  TO    Speed  40  Knots 

Det.  range  750  m  Lobe  width        20° 

Sweep  angle      30 
Turn    rate  187s 


O       0°   Attack   angle 

a     60° 

A  120° 


One  detection  — 
uccessive 
etections 


1500 


3000 


5000  7000  m 

Firing    range 


Figure  21. b.  -  EFFECT  CF  FIRING  RANGE 
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Torpedo    Parameters.  Tactical    Situation 

TO    Speed  32       40    Knots  TA    Speed      18  Knots 

Lobe    width        20      20    °  Det  range     750  m 

Sweep    angle       40       30    ° 
Turn   rate  15       18    °/s 


n 
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<d 
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One   detection 


o  1500  m  Firing   range 

o  3000m 

a  5000m 

A  7000m 


0  30°        60°        90°        120°      150°       180' 


Attack    angle 


Figure  22  - 


COMPARISION  OF  TWO  TORPEDOES  WITH  CHANGE  IN 
FIRING  RANGE 
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Torpedo     Parameters 


TO    Speed 

32 

40  Knots 

Lobe   width 

20 

20  ° 

Sweep   angle 

40 

30  ° 

Turn    rate 

15 

18°/s 

Tactical    Situation. 

TA  Speed      18  Knots 
Det.  range      750  m 


Two  successive 

detections 


o    1500m  Firing    range 

o  3000m 

a  5000  m 

A  7000m 


0°  30°        60°       90°        120°       150°      180° 

Attack   ancle 


Figure    22. b.     -    COMPARISION    OF    TWO    TORPEDOES    WITH    CHANGE 

IN    FIRING    RANGE 
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Table  VII  -   VABIATIQN  IN  FIRING  RANGE 
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J.   EFFECT  OF  TARGET  SPEED 


Generally/  we  anticipated  a  degradation  in  MOE  as  the 
target  speed  was  increased.  And  overall,  this  was 
confirmed. 

The  simulations  were  carried  through  at  12,  18,  24,  30  knots 
target  speed. 

However,  fig.  23. a  and  b  shows  some  interesting  patterns 
regarding  optimal  attack  angle  for  different  target  speeds. 
For  a  32  knots  torpedo,  at  60  degrees  attack  angle,  the 
torpedo  is  equally  good  for  any  type  of  target  speed  fcr  one 
detection  only.  For  two  successive  detections,  the  torpedo 
is  equally  good  between  30  and  90  degrees  for  12  and  18 
knots  target.  A  24  knots  target  gives  a  consistently  lower 
MOE  over  the  whole  range  of  attack  angles,  and  the  2 
simulation  runs  with  a  30  knots  target  confirmed  that  trend 
for  the  32  knots  torpedo. 

We  may  form  the  conclusion  that  for  one  detection  only 
60  degrees  attack  angle  is  an  optimal  attack  angle  for  the 
range  of  target  speeds.  For  two  successive  detections,  30 
to  90  degrees  attack  angle  gives  equally  good  MOE  between 
target  speed  of  12  and  18  knots. 

One  interesting  point  is   that   it   seems   that   if   the 
target   speed   is   less  than   0.4   of  the  torpedo  speed  the 
optimal  attack  angle  shifts  forward  to  0  degree. 
This  also  applies  to  two  successive  detections. 
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Tactical    Situation. 
Range         3000  m 
Det.  range    750  m 


Torpedo    Parameters: 

TO    Speed  32   Knots 

Lobe  width  20  ° 
Sweep  angle  40  ° 
Turn    rate  15°/s 


0°  30°  60° 

Attack  angle 

Figure  23  -   EFFECT  OF  TARGET  SPEED 
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Tactical    Situation: 

Range  3000   m 

Det.  range      750  m 


Torpedo     Parameters: 
TO    Speed         40  Knots 
Lobe  width      20° 
Sweep   angle    30° 
Turn   rate  18°/s 


One   detection 
Two  successive 
detections 


150"         180° 
Attack    angle 


Figure  23. b.  -  EFFECT  OF  TARGET  SPEED 
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For  a  40  knot  torpedo,  in  addition  to  the  point  of 
optimal  attack  angle  at  0  degree  for  slow  target  speeds,  we 
also  experienced  a  relatively  low  MOE  for  slow  targets  in 
the  range  45  to  105  degrees  attack  angle,  compared  to  fast 
targets.  But  as  a  compensation,  MOE  is  increased  for  small 
attack  angles  and  the  astern  attack  angle  compared  to  fast 
target.  Obviously,  some  type  of  a  breaking  point  is 
experienced  for  target  speed  of  . 4  or  less  of  torpedo  speed. 

Why  a  slow  target  produces  this  increase  in  MOE  in  the 
two  extreme  cases  (ahead  and  astern)  may  be  explained  by  the 
balance  between  time  to  reach  detection  range  and  the  total 
relative  speed.  It  is,  however,  more  difficult  to  give  any 
explanation  of  why  a  slow  target  should  produce  a  lower  MOE 
for  some  attack  angles  than  a  faster  target  does.  One  would 
have  anticipated  an  increase  in  MOE  over  the  whole  range  of 
attack  angles  for  a  slow  target. 
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Torpdo     Parameters:  Tactical    Situation: 

TO    Speed         32     40Knots        Range  3000m 

Lobe  widtn      20     20°  Det.  range      750m 

Sweep   angle    40     30° 
Turn   rate  15      18°/s 


One 
detection 


.8 


7 


>  .6 


n 
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a    60° 

A  120° 
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TA    Speed 


Figure  24  -   COMPARISION  OF  TWO  TORPEDOES  WITH  CHANGE  IN 

TARGET  SPEED 
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VII.   TACTICAL  ANALYSIS 


In  addition  to  the  detailed  parametric  analysis,  which 
has  been  shown  previously,  we  also  could  expand  the  analysis 
to  cover  more  tactical  related  problems.  If  we  assume  a 
given  target  speed,  we  could  construct  detection  probability 
charts  as  shown  in  Fig.  25.  a.  and  b. 
This  analysis  would  then  naturally  fall  into  two  areas: 

direct  comparison  of  two  or  more  different  types 

of  torpedoes. 

effect  of  tactical  situation  on  the  detection 

probability. 

The  two  charts  (Figs.  25. a.  and  b.)  were  formed  by 
running  simulation  runs  for  differnt  tactical  situations 
(range  and  attack  angle) ,  and  then  fitting  constant 
detection  probability  curves  through  the  data  points. 
The  use  of  these  types  of  charts  falls  into  two  areas: 
Evaluate  different  torpedo  types  for  different  tactical 
situations;  essentially,  which  torpedo  is  best.  Or  for  a 
given  tactical  situation,  how  could  the  situation  be 
improved,  and  what  options  exist. 

The  first  type  of  use  applies  mainly  to  operational 
planning;  operational  requirement  in  the  design  phase  of  a 
torpedo  and  procurement.  By  laying  one  chart  atop  of  the 
other;  we  get  a  visual  picture  of  how  much  is  improved  when 
using  a  'better'  torpedo,  and  for  which  tactical  situation. 
The  shaded  area  in  Fig.  25. b.  shows  how  many  more  tactical 
situations  have  been  covered  when  going  from  a  32  knots 
torpedo  to  a  40  knots  torpedo  for  0.25  detection 
probability. 
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The  other  type  of  use  of  the  charts  is  tactical.  When  a 
firing  unit  decides  to  attack,  and  finds  itself  in  a  given 
tactical  situation,  the  question  is:   What  to  do  ? 

For  given  target  speed  and  own  max  speed,  the  charts 
make  it  possible  in  a  simple  way  to  decide  where  to  go  and 
what  course  to  keep.  Also  from  the  charts,  one  can  decide 
where  on  the  relative  course  is  the  optimal  firing  position. 
For  a  submarine  attacking  a  zig-zagging  target,  the 
Commanding  Officer  can  better  make  his  evaluation  of  when  to 
fire,  as  the  attack  angle  and  the  distance  are  continuously 
changing.  He  can  see  what  improvement  to  expect  when  the 
target  will  change  course  next  time.  An  example  of  a 
tactical  situation  and  the  course  of  action  to  follow  are 
given  in  Fig.  25. a. 

These  points  also  bring  up  the  question  of  what  to 
improve  in  the  operational  picture;  the  firing  unit's 
ability  to  acheive  a  good  firing  position  or  the  torpedo's 
ability  to  detect  target  from  non-optimal  situations. 
In  this  discussion,  the  guided  torpedo  has  to  be  brought 
into  the  picture.  The  effectiveness  of  guidance  has  not 
been  adressed  at  all  in  this  study,  basically  because  that 
would  have  significantly  expanded  the  scope  of  the  study,  as 
well  as  bringing  in  the  whole  problem  of  fire  control 
equipment,  its  effectiveness  and  its  reliability. 
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VIII.   CONCLUSIONS 


The  study  was  carried  out  in  order  to  investigate  the 
detection  process  of  an  active  sonar  homing  torpedo  used 
against  surface  ships. 

Specifically,  we  wanted  to  study  the  effect  of  changes 
in  torpedo  parameters  such  as  torpedo  speed,  turn  rate, 
sweep  angle  and  detection  range,  as  well  as  changes  in  the 
tactical  situation  such  as  target  speed,  firing  range  and 
attack  angle. 

In  an  attempt  to  gain  insight  into  the  complexity  of  a 
homing  torpedo,  the  described  model  was  built  and  the 
simulations  done  as  previously  shown. 

In  designing  a  homing  torpedo  and  evaluating  torpedo 
tactics  the  detection  probability  is  an  essential  part  of 
the  total  effectiveness  of  the  torpedo. 

To  be  able  to  hit  the  target,  the  torpedo  has  first  to 
detect  it,  which  justifies  why  we  started  out  with  analyzing 
the  detection  process. 

Also,  as  part  of  this  analysis  we  investigated  certain 
aspect  of  the  next  step  in  the  operational  process; 
acguisition. 

It  is  not  difficult  to  visualize  tests  which  may  be  used  in 
order  to  recognize  an  echo  as  a  detection  and  subsequently  a 
target  to  attack.  Some  of  these  tests  may  be  doppler, 
successive  detections,  detections  within  a  given  range, 
2-of-3  detections,  size  of  echo,  length  of  echo  etc. 
The   problem   of   false  echo,  however,  was  net  approached  in 
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this  study.  That  would  have  to  be  the  next  to  consider  in 
relation  to  reducing  the  number  of  successive  detections  in 
order  to  acquire  a  target. 

In  order  to  allow  for  the  errors  in  tracking  of  the 
target  before  firing  and  also  small  maneuvering  of  the 
target  after  firing,  we  introduced  errors  in  the  target 
speed  and  course  when  calculation  of  the  torpedo  main  firing 
course  was  done. 

During  runs  the  torpedo  was  unguided,  and  did  not  react 
on  any  detection;  i.e.  it  did  not  attack  the  target.  For 
sonar  condition,  isovelocity  was  assumed  and  no  surface 
effect  was  built  into  the  model. 

The  result  gave  certain  insight  into  the  complexity  of 
the  detection  process,  stressing  the  importance  of  a  good 
tactical  firing  position,  and  of  a  high  speed  torpedo  with 
long  detection  range. 

However,  the  data  also  showed  the  relationship  between 
detection  range,  lobe  width  and  turn  rate,  as  well  as 
weighting  the  sweep  angle  in  relation  to  torpedo  speed.  It 
can  also  be  concluded  from  the  results  that  changes  in 
torpedo  parameters  as  turn  rate  and  sweep  angle,  which  may 
be  inexpensive  modifications,  will  not  give  a  significant 
improvement. 

Generally,  the  overall  important  factor  was  the  time  the 
torpedo  used  in  order  to  travel  within  the  detection  range 
of  the  target.  This  was  due  primarily  to  the  error 
generated  in  the  target  data,  obviously  the  actual  value  of 
the  result  is  sensitive  to  these  assumptions. 
However,  the  understanding  and  insight  in  the  detection 
process  achieved  by  simulation  should  not  be  reduced  by 
other  assumptions  with  regard  to  error  in  target  data. 
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With  regard  to  the  analysis,  the  result  has  shown  a 
consistent  and  general  trend  that  if  we  are  able  to  require 
only  one  detection  for  acquiring  a  target,  the  detection 
probability  is  significantly  higher  than  if  more  than  one 
detection  is  required.  And  what  is  more  important,  the 
potential  for  improving/optimizing  a  homing  torpedo  is  also 
significantly  higher  for  one  detection  only. 

This  implies  a  large  payoff  for   other   methods   of   keeping 
down  the  probability  of  false  detections. 

Secondly,  a  high  speed  torpedo  has  shown  a  general 
superiority  in  MOE.  This  was  specially  obvious  in  attack 
angles  greater  than  90  degrees,  which  tends  to  make  a  high 
speed  torpedo  more  of  an  all-round/reliable  torpedo  with 
regard  to  tactical  situations. 

Thirdly,  except  for  changes  in  attack  angle  and  firing 
range,  the  detection  range  seems  to  influence  the  MOE 
strongly. 

These  three  remarks  all  point  towards  an  improvement  in 
the  sonar-/filtering-area  as  the  most  promising  area  in 
which  to  carry  out  research  and  invest  effort. 

This  study  has  also  pointed  out  the  advantage  of  high 
torpedo  speed  and  firing  at  short  ranges.  There  exists 
therefore  considerable  argument  for  a  short  range,  high 
speed  torpedo,  given  that  one  is  able  to  position  the 
torpedo  at  a  short  firing  range;  i.e.  a  small,  simple 
torpedo. 

Basically,  there  are  two  schools  of  thought; 

a  highly  sophisticated  torpedo;  long  range,  guidance, 

expensive,  but  close  to  the  one  shot-one  hit  idea, 
a    simple,    high   speed    torpedo;     short    range, 
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non-guidance,  inexpensive,  and  requiring  either  a 
firing  unit  which  can  get  into  an  optimal  firing 
position  or  a  larger  number  of  shots  to  acheive  hit. 

The  result  may  be  useful  in  giving  example  of  how 
tactical  guidelines  can  be  evaluated  by  the  simulation 
approach.  3ut  more  significant  is  pointing  out  the 
importance  of  torpedo  capability  and  the  tactical  situation. 
Obviously,  we  have  to  look,  on  the  whole  torpedo  system, 
including  the  firing  unit.  Investment  in  resources  and 
effort  should  not  necessarily  be  spent  only  on  the  torpedo 
in  order  to  increase  its  effectiveness,  but  may  be  spent  on 
the  firing  unit  as  well  in  order  to  make  the  unit  able  to 
reach  a  tetter  firing  position. 

A  follow-on  of  this  study  may  be  to  investigate  the 
attack,  process  of  the  torpedo,  including  the  acquisition- 
and  hit-problem. 

Then  the  question  of  guidance  during  torpedo  run  should  be 
analyzed  in  order  to  better  evaluate  the  problem  of  choosing 
between  a  few  sophisticated,  expensive,  guided  torpedo 
system  or  many  simple,  inexpensive,  nonguided  torpedo 
systems. 
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APPENDIX    B 


FLOW    CHART     FOR    SIMULATION     PROGRAM 


A    TOREEDC    SIflULATICN. 


aAIN     EHOGEA3 


PAG2     1 


I 


A  ICHPECO  ilJULAIION.       HAIN  PHOGRAa. 

A  TORPEDO  SiaOLAIICN. 

SiaULAIING  AN  AC1IVE  HONING  TORPEDO  CUBING  SEARCH. 

THE  PROGRAM  IS  HUM  IN  0.3  SEC  STEPS. 

CGflaON  ISEED2,  TTiaE,  TO,  IA ,  TSAI2,  SANGE,  AiFA,  LAaBC,  TADEC, 
BEAR,  RAE,  TAC,  CCCR,  3EVSP,  ENG,  EM,  PH2,  aCCURS,  TCOORS, 
BXT,  aXH,  IK,  IDTI«2,  ITia2,  XT,  YT,  XTAR,  ITAR,  TDIST,  .1DIST, 
TUENTO,  ISTVA1,  PHI,  RMAX  (5  ,  6) , TRANGE ,  OIST.IPRINT 

CCaaON/DAIA/EA,OPLOB 

CGaaOM/TAHG£T/TACaG,IAai ,BNGaOD ,DA1,COOR,CE (10)  ,SE(15)  ,JHUN, 
If LAG,2A2 

DIMENSION  Ba(5,5),  VAR(5),  DET(150,5),  DETB  j 1  50,  5)  ,  STD  (5)  , 
ASE  EC  (15  0,5)  ,DERB(  150,5)  ,CLOS3{  150,5)  ,<CON  (150,  5f 

REAL  LAflEE,  aCOURS,  HIS,  axa,  aDIST,  LAflBDG 

INTEGER  50NOUI 

SETTING  Of  CONSTANTS  (STEP  1) 


CAI1  OVFLCSi. 
j_ 

EHI   =3.iai59265H 
EH2   =2.*P9I 
SAD   =EH2/360. 
ISEED1  =262776 
ISEED2  =961695 


I 

CCOREC   -   BAX  2RBOB  IN  TABGiT  COUESE  ESTIBATE 

_i_ 

I  CCOREC  =15.  I 

CCOH  =CCOBEC*aAD 


CSPEED   -   SI  DE7  IN  TABGET  SPEED  ESIIMATE 
I 


CSPEED  =  3. 
C£VSE=CSE£E£*0.5 
IETiaE  =1 


I 


SET  NOaEES  OF  ITERATIONS 


I 


I  IRUN  =150 


(CONTINUED  ON  EAGE    2) 
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A  TOEFEDO  SIHULATION.      SAIN  PHQGflAa.  PAGE 

I 


SET  PHINT  OUT  HGDE 

_i_ 
1  IPBIBT  =T 

I 


SET  LCEE  0?F  IOREEEO  CENIEB  3EA2ING 

I  CPLO5=0.  I 

I 

SET    TAE1ES    TC    ZEEO     (STEP    2) 

JL 

♦  DO  + 

♦  +  +  ++>+-»•■»+  11  ♦ 

♦  +    1=1, 5  + 

+  i 

♦ 

+  +--- + 

♦  +    DC  + 

•f +  +  +  +  ♦-••■••++  12  ♦ 

♦  ♦    J=1,5  + 

♦  I 

♦  \_  

♦  |     5M(I,J)     =0.  I 

♦  I 

12  ♦-♦♦+-t-+++CCNIINUE 
+ 

11    -f-f++++++CONTIN0E 

I 

+- * 

♦  DO  + 

+  •*  +  -♦+♦+  +  +  +  13  + 

♦  ♦    J=1 , 5  + 

♦  + + 

:  I 

♦  + + 

♦  +  co  ♦ 

+  •«•+  +  ♦■♦•++-«■+  14  + 

+  -*■    I=1,IEUN  + 

♦  I 

+  __[_ 

+  I     KCM(I,J)    =0  I 

♦  I 
14    +♦+++++ +CC STINOE 

13  +  +  ++*  +  +  fCCJiTINOE 

I     EA2      =0.  I 

(CONTINUED    ON    eAg£         3) 
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4    TOBP.EDC    SIMULATION. 


MAIM    PROGBAH. 


PAGE 


COHPOTE    TARGET    EBSOBS    AND    STOB2 


*    DO 

♦  ■♦■♦■♦■♦*  +  ♦+  +  20 

♦  +1=1, 10 

♦ 

♦ 
•f 


I    CE(I)=-  (15.*BAD+CCOH/10.)  +  (CCOB 
10.) *!*? 

I 


20  ++*+++++CCSIINU2 


100  ♦    DO 

♦  ♦■«■♦+••■+  +  -•■  +  900 

♦ 
•f 
+ 


+  JBUN=1,ISaN   ♦ 


SET  SON  COUNTERS  (STEP  3) 


IELAG=1;  TOC  IOH  TCEPEDO  5PEZD 
I 


IPIAG=0 
ICCNT=0 
JCONT=0 
ITI3E=0 
1L  =C 
IK  =0 
CIST  =0. 

axT  =0. 

3X3   =0. 

SUNCUT  =0 

I 


(CONTINUED  ON  PAGE 


*) 
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A  TORPEDO  SIMULATION.  BAIN  PROGRAM.                               PAGE 

♦ 

+ 

♦  REAL  IN  SETTING (TORP  AND  TACTICAL) 

♦  FIRST  HON  ?  (STEP  4) 

:  .  J  * . 

♦  .  *     IF  *  . 

♦  *  JHON.GI.1      * 

♦  *  .  .  * 

♦  *  .  .  *          T     |  160   | 

:  ']' 

♦  I 

♦  I 

♦  CALL  PARMET 

+ 
+ 

♦  ***WRIIE  (6,228) CFLOB 

228  *  FORMAT  (1X,/.  1X,*SCN1-S  MAIN  LOSE  OFF-SET  FROM  CENTER  ', 

+  ■  SEARING',  F6. 2,'  TIMES  DEFLECTION  A3GLS*,/) 

♦  _j__  

+  I  10151=10/2.                             I 

♦  j  MDIST=TA/2.  I 

♦  I  IURNTC  =TEATS/2. 

♦  I  INTVAL  =IFIX((TTi:iE/C5)  >0.5)  | 

♦  i 

♦  INTVAL  GI7ES  NOHEER  OF  IIMESTSPS  FOR  EACH  TRANSMISSION 

+ 

♦  i 

♦  .  *     IF  *  . 

♦  *     IPRINI.EC.O     * 

♦  *  .  .  * 


.  *         T     |  160   | 

♦  *  .  *  

+  I 

♦  F 

:  I 

♦ 

♦  PRINT  CF  HEALING 

•f 

151    ♦  ***HRITE (6, 220) 

+ 

230    +  FORMAT  MX,  // ,  1X,  'RUN'     4X,  '  EST    OF    TARGET'  , 7X  ,' TORP       TORE    M' 

f  8X,'TCRP    CCCSC',5X, 'TARGET    COORD            RUN* , 31 , ' TOR?' , / , 

♦  1X,'NO  CG0R5E    SPEED       RANGE  DA • , 4 X ,  'COURS E '  ,  1  OX  ,  '  X  '  ,  6X , 

♦  ,X'/7X,,X,,6Xf,I,,6X,' STOP'  ,2X,  '  RON' ) 

(CONTINUED  ON  PAGE    5) 
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A  TOHPEEO  SIMULATION. 

♦ 

♦ 
+ 


MAIN    PB0GBAM. 


PAGE 


♦ 

160    ♦ 

* 
* 
+ 

+ 
♦ 

+ 

+ 
+ 

* 
♦ 

♦ 
+ 
+ 

♦ 

+ 

+ 

+ 
♦ 
•f 
+ 

198    ♦ 

+ 

♦ 

♦ 

+• 

152    ♦ 


CALCULATE    TOHPEEO    STABT    JOSITION     (STEP    5) 


XTAH    =1500C. 
YTAfi    =1500C. 
CPA1    =EANGS 
£NG       =  TAC*EEA.B 


I 


*  IE 

ENG.G1.PH2 


|     ENG=3NG-PH2 


I 
.     *         *    . 

*  IE 

3NG.L1.C. 


I     3N=PH2*BNG 


XI         =  XTAB  +  HANGE*SXB(BNG) 

21  =YIAB+RANGE*COS <BNG 


*  IF  * 

IPHIM.EC 1 


I     152       | 


***WRII2(6,  IS  6)  JBON 

?OSflAT(U,/,6X,'RaN    30HB2B 


:'  ,IU) 


CALCULATE    IOSEEEO    EEFLECTICN    ANGLE    AND    EI3ING    SITUATION 

CALL    PIBIHG 
(CONTINUED    ON    PAGE         6) 
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SAIN    PROGBAtf. 


A  TORPEDO  SlflOLAIION. 

+ 

♦  SEI  D2TECTI0NTAS1ZS  TO  ZEBC  (ST2P  6) 

+ 

♦  _J 

♦  ♦  DC  ♦ 

♦  ♦♦■♦♦♦♦♦♦  +  10  * 

♦  +    1=1,5  ♦ 
+                           +. _  + 

♦  I 

*  — L. 

♦  ♦    DO  * 

♦  ♦♦+-••♦•♦•+++  9  ♦ 
+                        ♦    J=1,6                      ♦ 

♦ 
+ 
♦ 


PAGE 


I     BHAX(I,J)     =C, 


9    ♦  ♦♦♦♦"•"►♦CCNTINOE 
10    ++++>+++CCSTINUE 


♦ 
•f 
♦ 

+ 
•f 
+ 
+ 

♦ 
+ 

+ 
■f 

+ 
•f 
+ 

:31     ♦ 

+ 
♦ 

+ 


TEST    TOBPEDOSCN    CDT     (STEP    7) 
I 


I     THUS    =TEANGZ/TO 

i    CLOSP=TA*CCS (BEAB)  *T0*COS (DA) 

|     HBUN    =  aA!JG£/CLOSP 


-     *  Is  * 

*         IB0N.G1.HB0B 


*    .    * 

i 


*— -- 


|     499       | 


***H5iis  (6,  ::i) 

FCBflAT(1X,/, 1X,'TABG2T  OUTSIDE  TOBPEDO  BANGE') 
Z2BC  TAEI2S 


♦    DC  * 


♦++++♦++++  23 

♦  ♦    J=1,5 


i     CITB<JHUN,0)     =0. 


(CONTINUED    Oil    PAGE         7) 
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A    TORPEEO    SICOIAIICS.  SAIN    PROGRAM.  PAGE 


,i)=o. 


♦  DEI  (JHON 

♦ 

♦  - — 

♦  I     flSPEC(J2US,J)    =0.                                                           | 
+  1     CI0S3(JSUN,J)     =0. 

+  I     EEBE(JBON,J)    =0.                        •                                     I 

♦  I 

23    ♦  ♦♦♦♦•♦■♦  +  CCSTIUUE 

♦  I 

+  

♦  I     9C0     | 

+ 

♦  t 

♦  1 

499  ♦  .    *         *    . 

♦  .    «            IE  *    . 

♦  *            IELAG.2C..1  * 

♦  *    *  .     *                         T             |     990       | 

♦  *    .    *  

♦  I 

♦  E 

♦ 

♦  I 

500  +  .    *         *    . 

♦  .    *            IF  *    . 

♦  *  aXl.GI.TSd.NGZ  * 

+  *    .                    .    *                         T            |     aOHOUT=1 

♦  ••»    .    *  

♦  I  I 

♦  f  I 

J  I 

♦  .    *       *   . 

♦  .    *         ii  *   . 

♦  *  RUNO01.EC..1  * 

♦  *    .  .     * 

+  *    .                     .    *                         T            |     990       | 

+  *    .    *                                                 

:  i 

:  I 

+ 

♦  CALCULATE  SEH  POSITIONS 

+ 

+  CAIL    PCSIS 

♦  CHECK    IE    TARGET    IS    DETECTED 

<CCNTI.VIGED    Oil    eIgE         8) 
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A  TOBPEDO  SIBULATICH.      3AIN  EHOGflAfl.  PAG2 

♦  CAI1  DETECT 

+ 

♦  CHECK  CPA (CLOSES!  POIHI  DP  APPROACH)   (STEP  3) 

♦ 

♦  _j_ 

♦  I    Ii         =11  +  1  I 

♦ 
+• 

♦ 

♦  *    .  .    *  T  |     500       | 

♦ 

+ 
•f 
+ 

155    ♦ 

♦  *    .  .    *  T  |     RUNOOT=1 

♦  i  i 

♦  __  .j. 

♦  i 

♦  |     CPA1    =CEA 

:  T 

♦  

♦  I     500     | 
+                                  

♦ 

♦ 
+ 

♦  GENEBA1E  STATISTICS  (STEP  9) 

990  +        CCNTINOE 

♦  ♦    CO  + 

♦  +  +  +  >+-♦•+■♦•+  910  ♦ 

♦  ♦    IK1=1,5  * 

♦  I 


* 
a 

.J 

.  *        If 

IL.LI. 

*    „ 

*    . 

,20 

* 

* 

* 

I 
I 

I 

J 

IL 
CPA 

i 

=  0 

=  CIST 

I 

——. 

* 
* 

* 

CPA1- 

I 

*    ! 

II 
•CPA. 

,L£. 

0. 

(CONTINUED  ON  PAGE    9) 
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A  TOBPEDO  SIflOLATICN.      MAIN  FROG  BAH.  PAGE 

♦ 
•f 
♦ 

+         KON<    )     =    DETECTICN/NO    DETECTION 

♦  .   *       *  . 

♦  .   *         IE  *  . 

♦  *    HfiAX  (IKI,1)  .GT.1.     * 

♦  '    *    .  .    *  T  |     KON (JBUN,IKL)=1 

♦  *    .    *  

:  i I 

+ 

♦  DET  <  )  =  DISTANCE  TO  TARGET  AT  DETECTION 

+ 

♦        \_ 

♦  |  LET  (JH0K,IKI)  =  E2AX  (IKL, 1)  I 

♦  I 

+ 

+    DETE(  )  =  3EAEING  TO  TAEGET  AT  DETECTION 

*> 

♦         _l_ 

♦  I  EET3(JR0N,IKL)  =  BMAX  (IKL,  2)  | 

♦        — — 

+  I 

♦  ASPEC  (  )  =  TAEGET  ASPECT  AT  DETECTION 

♦ 

♦  ASPEC{JR0S,IKI)=HtfAX.<IXL,5) 

+    I 

♦ 

♦  CICS3<  )  =  5EJBING  TO  CLOSEST  PABT  OF  TARGET 

+ 

♦  CLCS3  (JHQH, IKL) =2aAX(IKL, 4) 

+ 
■f 
♦ 
•f 

♦  DE2E  (  )  =  BEL  EEABING  FROM  SAIN  TORP  COUBSE  TO  TABGET 

♦  l_ 

♦  |  CEBE  (JSUN,IKL)  =EMAX(IKL,6)  | 

♦  I 

910  ++++++++CCSTIS0E 

♦  I 
{CCNTINOED  ON  PAGE   10) 
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A  T05PED0  SIflULAIICN. 

+ 

+ 

+  * 

+• 
♦ 
♦ 
■* 
•f 
■► 
♦ 

+ 


MAIN    EBQGEAil. 


PAGE       10 


.    *         *    . 

*  If  * 

IFLAG.EQ. 1 


I    991       J 


+ 
♦ 

+ 

+ 
+■ 

+ 
232    * 

■¥ 
♦ 

■f 
♦ 

♦ 
•f 

+ 

+ 
+ 
♦ 

195    ♦ 
200    ♦ 


PRINT    CUT     (SUP     TO) 

|  COUR  =SCGOFS/RAD 
!  IC  =:CCD55/RAD 
|     II         =IlIHE/2 


*  IE 

IPBIiiI.EC.O 


I     195       \ 


***«RITE  (6,  222)  JRON  ,  TAC^IG  ,IAM1  ,  3NGaGET,DA1  ,COOa,XT,YT,  XTAa,TTAa, 

FCEiiAT  (1X«/.  1X,I3,3X,F5.  1,2X,F<*.  1,  1X,P6.0.  2X,F5.  1,32, 
F5.  1,7x\  ft.  0,  1X,Fd.O.,  2X,F6.0,1X,F6.0,  3  2,14  ,2X,F6.0) 

I 

.     *         *    . 

.    *  IE  *    . 

*  EBAX  (1,1).  2C..0.        * 

*  .  .     * 


991 


I     9C0     | 


***«HIIE  (6,200)  IT 


FCEMAT  (12,'aUH    STOPPED    AFTER     ',15,'     SECONDS',//, 
1X,'3UN    DATA    AS    FOLLOWS    AT    END    OF    HON') 


(CCNTINOED    ON    PAGE       11) 
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A  TOBPEDO  SIMULATION.  3AIN  PROGRAM.                                PAGE  "  11 

♦  I 

■* 

♦  ***WHIIE(6,2C2) SXT,DIST 

202    ♦  PORflAT(lX, 'TOTAL    TOflP    HUN    •  ,  F9.  1,/,  1  X,  •  DIST    TO    IABGZT     «,F9.1) 
+ 

♦  ***WBIIE  (6, 2C4)  2T,YT,XT-A2,YTAR,C0Ua,TC 
♦ 

204    ♦  FOBMAT  MX.'TOBP    X-COQRD    ■ , 2X  ,P9 . 1 ,4X, ' TOBP    Y-COORD     ',2X,F9.1,/, 

♦  1X,'TAflGEl    X-CCORD    • , F9. 1 ,4 X , ■ TABG2T    Y-COORD    «,F9.1,/, 

♦  1X,'TORP    flAIN    COURSE    •  ,F9 . 3, 6X, ' TORP    COURSE    ',1Z,F9.3) 

:  I 

♦  .   *       *   . 

+  .     *            IF               *    . 

♦  *       EMAX  (1,  1)  .F.Q.O.       * 

♦  *    .  .     « 

♦  *  .  .   *                  r         i   991     i 

+  *   .   *                                   

:  t 

:  ! 

♦  «**«EIT2  (6,206) 

+  ***WBITE  (6,2C8)   (I,  (RMAX,(I,J)  ,J=1,5)  ,1=1,5) 
+ 

206    +  FCEflAT  (1X#/, 1X,' MAXIMUM    DETECTION    RANGES    AND    BEARINGS', 

♦  /,  1X, 'SUCCESSIVE'  , 4X,' SAX       DBT1 ,9X,»DEX    3EABING' , 


+ 


/,  1X, 'SUCCESSIVE'  , 4X,' MAX       DEI  •  ,9X,  »  DF.T    3EABINI 
5X,'MAX       CET'  , 11X, 'SET    3EARI NG ', 4X ,' TARGET • ,/, 


♦  1X,'D2T  NC.',4X, ' aANGE    -    CENTER', 4X, 
+  'CENTER*  ,SX,  'RANGE    -    CLOSEST'  ,4X  ,'  CLOSEST  '  , 

♦  3X, 'ASPECT') 
+ 

208    f  FOBMAT  (32  ,14 ,7X, F 1 S . 1 , 7X , F3 . 2  ,7X,F10. 1 ,  9  X, 

♦  F8.2,7Xf F6.2) 

♦  I 

♦  |     900     | 

+ 
+ 
■f 

991    i-  ***WEITE  (6,220) 

♦  ALL    BOS    COMPLETED    ?     (STEP     11) 
SOO    *>>*++++CCNTINUE 


CALCULAIE    SUMMASY    RESULT     (STEP    12) 


♦    DO  f 

♦  ♦■♦♦♦♦++-f  ■»  27  ♦ 

♦  ♦    KH=1  ,5  ♦ 
+                           + , _ — + 

♦  I 

♦  i 
(CONTINUED  ON  PAGE   12) 
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A  TOBPECO  SiaOLATION.       MAIN  PROGBAil.  ?AG2   12 

I 

+  ♦—— ♦ 

♦  +    DO  ♦ 

♦  ♦♦-♦•♦♦♦+++  28  + 
+                       ♦    L8»1,IB0H           ♦ 

+  •*• + 

:       I 

♦  I     Ba(KB,1)     =  RB  (KB, 1)  *?LOAT  (KGN(LR,KB) )        | 

♦  £a(KB,2)    =  Efl  (KB,2)*DET  (LB,KH) 

♦  |     Ea(KR,3)     =FB(K2.3    ♦  ASS (ASPEC  (LR , KB) )        | 

♦  |     Ha(K2,4)     afia <KB,4    +DET3 (LB,K2)  | 

♦  I     FH(Ka,5)    =  EB  (KB,  5)  «-DEa3(LB, KB)  | 

♦  I 

28    ♦♦♦♦++++CCKTIHOE 
27    >+  ++>+-f+CCSTINl/E 


♦ 

♦  I     EEL      =SE(K?,1) 

+ 

+ 

+  .     *         *    . 

♦  .     *  IE  *    . 

♦  *  DEL.  12.1.  * 

♦  *    .  .    * 

+                            *    .                     .    *                         I             |     DEL=1. 
«.  *    „     *  

:  i  I 


♦ 

♦  t  5B(KK,2)     =  EB  (KR,  2)/DEL 
+  I  RB  (KK, 2)     =SB    KK,3    /DEL 

♦  I  £3    KK,4)     =EB    KK,4    /DEL 

♦  l  5H(KK,5)     =  nfl  (KK,  5|  /DEL 

♦  I 
920    +-*+++>  +  +  CCSTINUE 


(CCNTINOEE    ON    SAGE       13) 
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A  TOHPEDO  SIMULATION. 


HAIN    EHOGEAil. 


?AG2       13 


* — — — + 

♦    DO  + 

■*■*■*■  +  ■*+++  +  +  25  ♦ 

♦  ♦    XX=1,£  ♦ 

:  T 

♦  -— _- — _ 

♦  I     VAB(KK)     =0. 

+  I 

♦ 

♦  ♦    DO  * 
+++++++++♦          170  + 


♦    I*1,IB0H 


.    *         *    . 
.    *  I?  * 

*      RCS(I,BR)  .fC.O 


+ 

+ 
+ 
+ 

+ 
+ 

+ 
♦ 

170    ♦♦++++++COMIIHOE 

+ 
25    >  +  +  *+  +  ++CON'IINaE 


-— ♦ 


I     WO       | 


I     VAH(KK)     =7AB(KK)  ♦(  (£11  (I,KK) -fifl  (KK,  2)  ) 

I  **2) 


+    DO 
♦  4-*  +  -f-f+++ +  26 


+    KB=1,5 


|     EEL       *3£(KS,1) 


I 


*  12  * 

DEL. IE. 2. 


I    D2L=2, 


|     STD(KH)     =  SC,3T  (VAa  (KS  )/ (DEL-1 . )  ) 
J    Fa(0,1)     =?S  (XE,  1) /FLOAT  (IHOS) 


(CONTINUED  ON  ilGl       14) 
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A  TORPEDO  SICULA1I0N.      BAIN  PEQGRAH.  PAGE   14 


F8IJJT  SOHHABY  (SXEE  13) 

***WRITE(6,  1<7) JfiUN 

1S7  FCEaAT(1X,//,6X,'Sa;iHA£Y    0?    RESULT    AFTER'  ,  3X,IU  ,2X  ,  •  RU  NS  • ) 

*** WRITE  (6,  19  9)   (  (Hfl(I,vJ)  ,  J=1  ,2)  ,STD  (I)  ,  (RM  (I,  J)  ,J=3,5)  ,1=1,5) 

199  FCEMATHOX,  •  EECEAEILITY    OF    DETECTION' , 3X . 

■  AV<£RAGE',6X,rSTD    DEVIATION'  ,  7X  ,  '  A  VER  AGE  •  ,  7  X,  '  AVER  AGE  •  , 

72  'AVERAGE'  /  412 

•CET  RANGE'  ,bx''DET  RA NGE ', 7X, • TARGET  ASPECT'  ,4X ,' DET  BEARING', 

3X,'REI  EBABIS6' 


END 


220  ECEaAI(lX,»NC  DETECTION  MADE  DOHIMG  THIS  HUN') 

l_ 

I     CA2      =CA2/rICAT(IRUN)  I 

I 

***WHIIE(6,  TSO)  EA2 

190  F0BMAI(1X,/,  1X,' AVERAGE  DEFLECTION  ANGLE  :  •  ,52 ,  F3.  'i  ,/) 

**»HBITE(6,22U) 

234  FCBHiT(1X,/. 12  • DIST3IHUTI0N  OF  HUN  RESOLT  -  CENTER  OF  TARGET',/, 

6X,'0NE  SUCCESSIVE  DET ECTION  '  ,  1  OX,  'TWO  SUCCESSIVE  DETECTIONS', 
8X, 'THREE  SUCCESSIVE  DETECTIONS',/, 

2X,' BEAR'  ,22  'RANGE   ASPECT  BEAR  CLOS ' ,6X, ' EEAR' ,22 , 
'RANGE   ASPECT  EEAR  CLCS '  , 5 X  ,  •  BEAR • , 2X , 
'RANGE   ASESC1  3EAH  CLCS') 


***«H  IT  E  (6.236)  I  (DETB  (1,1), DET  (1,1)  ,  ASP  EC  (I,  1),CL0S3  (1,1)  , 
EElr  (1,2)  ,D2I  il,  2)  ,  ASPIC  (1 ,  2  )  ,CL0S3  (1 ,  2)  ,  DETB  (1 ,  3)  , 
DEI  (1,3)  ,  AS g EC  (1,3)  , CLCS 3(1,3) )  ,1=1, IRON) 

236  FOE  MAT  (3  (IX,  56.  1,  1X,F6.  1,  12,  ?6.  1,2X,F6.  1,5X)  ) 

999  STOP 
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I 
A  TOSPEEO  SIHOLATICN.     SUBROUTINE  PARHET. 

SOEEOUTINE  EiEaET 
BEAEING  IH  DATA  AND  PABAaETEBS 

COCMON  ISEED2,  TTI3E,  TO,  TA,  TBATE,  BANGE,  AL?A,  LAaSE,  TADEC, 
EEAB,  RA£,  TAC,  CCOR,  DEVSP,  BUG,    EN,  PH2 ,  SCCURS,  ICOUBS, 
MIT,  3X3,  IK,  IDTI2S,  ITiaS,  21,  YT.  ZTAR,  YTAR.  TDIST, 
aDIST,  TUENTO,  INTVAI.,  PHI,  EHAX(5,t>),  TRANG2,  Z>ISr,IPSINT 

BEil  LAHEE,  SCOUBS ,  i!2T,  3XM,  HDIST,  LAJ».BDG 

ISIEGEB  SUNCUT 


TEDEC  -  TECHNICAL  DETECTION  RANGE  (STEP  A1) 

LEVEL  CP  VARIATION:  37  5-750-1125-15  00  SETEBS 

TAEEC  -  IACTICAL  DETECTION  BANGE 
I 


IEDEC=750. 
IADEC=TEEEC 


ITIBE  -  TBANSIISSICN  INTEBVAL 

_l_ 

I  ITiaE=2.*IEEZC/150O.  | 

I 

TCKN  -  IOBP  SEESE  IN  SNOTS,  TO  -  TOBP  SPEED  IN  VSEC  (STEP  A2) 

LEVEL  CI    VABIATICN; 24-32-4Q  KNOTS 

__._    _J_  

I  TCKN  =40.  | 

|  TO    =TCKN/2  I 


TAKN  -  TARGET  SPEZE  IN  OOTS,  TA  -  TABGET  SPEED  IN  VSEC  (STEP  A3) 

LEVEL  Cr  VARIATION:  12-18-24-30  KNOTS 
I 


I  IAKN  =1£. 

I  TA    =TAKN/2 

I  . 
I 

(CCNTINQ2D  ON  FaGI    2) 


13U 
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I 

TACG  -  IABGZ1  CCCBSE  IS  DEGREE  (STEP  A4) 

r 


|  TACG  =270 

|  TAC   =IACG*BAD 


ALEAG  -  SWEEP  ANGLE  IN  D.SGS2E,  ALFA  -  S3EE?  ANGLE  IN  RADIANS 

(STEP  A5) 

LEVEL  Cf  VARIATION:  20-30-40  DZGHZES 
_j__      

|  ALEAG=30.  I 

ALPA  =ALEAG*2AD 


LABBDG  -  LOBE  WIETE  EACH  SIDE  0?  TORP  HEADING  (STEP  A6) 
IEVEL  CI    VARIATION;  10-20-30  DEGREES 


|  LAflEDG  =20. 

|  LAMED=LABB£G*RAD 


EELERG  -  aELATIY*  3EARING  FROM  TARGET  TO  TORP  IN  DZG22E     (SIE?  A7) 

LEVEL  CF  VARIJTICN:  0-30-60-90-120-180  DEGREES 
I 

|  BEL2BG  =-60.  I 

|  EEAR  =SEL£SG*RAD  | 


RANGE  -  DISTANCE  BETWEEN  TARGET  AND  TOBP  (STEP  A8) 
IEVEL  CI    VARIATION:  1500-3000-5000-7000  aETERS 

|  BANGE=3000.  I 

I 
(CCKTIN8ED  ON  PAGE    3) 
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I 

TBANGE  -  HA2  105E  FOH  IN  SETERS  (ST2P  A9) 

_       _j__   

|  TBANGE  =18CQ0.  | 

"1 

IEAIEG  -  TCRE  TUENEATZ  IS    DEGREE  PES  SZC  (STEP  A10) 
LEVEL  Of  VARIATION;  3-6-9-12-15-18-21  DZGRZ2/SZC 

I  IBATEG  =18.  I 

j  1RAT2=TBATEG*RAD  I 


CALCULATE  WIETH  CE  TACTICAL  SWEEP-LANE  (THZ0B2TICAL) 


1  SBNG  =TACEC*SIN(AI-?A+LAa3D)  *2-  | 

i 

CALCULATE  CCVEBA€Z  RATIO  (THEORETICAL) 
I 
I  CEAIIC  =1.- (TRAIE*ITIM2/(2. "LAMED) )     | 

I 
EfiINT  CF  SITCATICN  AT  START  0?  RON  {STEP  A12> 


.  *    *  . 

.  *     IE      *  . 

*  IP2IKT.EC..0     * 

*  .  .  * 

*  .        .  *         T     |  90 

*  .  *  


***WEIIZ (6,110) 

110  ECESAT  (11,//,  1X,  'TACTICAL  SITUATION  WHEN  FIRING', 6X, 

'TCRPEDO  EABAMETSS',/, 

22, 'RANGE   ATTACK   TARGET   T ARGZT ' ,6X, ' TEC . DET   T03P", 
3X, 'SWEEP   LC3E    TORN   SWEEP    COVERAGE',/, 
91, 'ANGLE    CCUHSE   5 f EED« , 7 X , « 3 ANGE     SPEED',2X, 
'ANGLE   WIDT6   RAT2   L ANZ' , 4 X , ' RATIO • ) 

(CONTINUED  ON  PAGE    4) 
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t 

***WRITE(6,  1  12) EANG2,BE1B£G,TACG,TAKN,TEDEC,T0KN,ALFAG, 
LAaECG,TBiTEG,SaNG,CSATIO 

112  FOB MAT  (1X.P6.0.22.F6.  1  ,3Xt ?5. 1 , 3X.F4. 1,61,16.  1,31,? 5. 1 ,31, 

FU.1,3XVf5.  1,22,?4.1f2X,7b.1,3X,F5.3) 

1*95   T 


90  **»»BITE(6,  100) 

100  FCEaATMX.'TACIICAL    SITUATION    WHEN    FI3ING ' ,/, 1X , 


FCEMAI  MX.'TACIICAL    SITUATION    WHEN    F: 
'FIBING    HANGE' ,31, • AITACK    ANGLE*     31, 
•TiEGET    CC0BSE',2Xr 'TARGET    SPIED*) 


♦»*WBIT2  (6,  102) RANGE, B2LE8G, TACG , IAKN 

102  FCEaAT(1I,4<2X,F6. 1,7X) ,/) 

***WBIT2  (6„  >C4) 

104  FCBBAT  (1 X,'TCBPEDC  PA5 AMETEBS • ,/, 1Xr 

•TECH.DET.2ANGE«,2X,.'  THANS  .1  NT.  V  AL  '  ,2X,'TOHP  SPEED'  ,  3X, 
'SSE2P  ANGLE') 

***«HIIE  (6,  1C6)  T2E2C,TTiaE,T0KN,ALFAG 

106  FCEMAT ( 1 X ,2 (21,21. 2,7X) ,2  (?6. 1 , 7X) ,/) 

***»fBIIE(6,  1C3)  LAaBDG,T2AI2G 

108  FOEMAIMX. 'LOBE    9IDTH « ,61, ' TU2N    BATS',/, 
3X,F6.1,  10X,F6. 1) 

***WBITE (6, 1C9) SBSG,C2AIIC 

109  FCEaATMX,'TiEC2ETICAL    WIDTH    OF    TACTICAL    S UEEP-LANE • , F9. 1 , 
/, 1X# 'TEEC2EIICAL   COVEBAGE    3AII0«,F9.4) 

95  2EIUBN 

END 
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A  TORPEDO  SIMULATION.    SUERCUTINE  FIRING. 


PAGE 


SUEROUTINE  FIRING 
CALCULATE  THE  TOEP  DEFLECTION  ANGLE,  MAIN  COURSE,  FIRING  CCURSS 

BASED  CN  ESTIMATE  OF  TARGET  DATA  (UNCERTAINTY) 


DIMENSION  0(2) 

COMMON  ISEEC2,  TTIME.  TO,  TA,  TRITE,  RANG2,  ALFA, LAMBD  .1 ACEC , 
BEAB,  SAD,  TAC,  CCCH,  D2VSP,  3NG,  EN,  ?H2,  MCCURS,  TCOURS, 
KXT,  MXM,  IK,  IDTIM2,  ITIM2,  XT,  IT,  XTAR,  YTA2,  TDIST, 
MDIST,  TUBNIC,  INTVAi,  PHI,  SMAX  (5,6)  ,  TRANGE,  DIST,I?fiINT 

CCMMON/DATA/EA,OFLOS 

CCEMON/TA5G2I/TACaG,IAM1 ,SNG MOD ,DA1 ,COUR, C2 (10)  ,SE(15)  ,JBUN, 
If LAG,DA2 

REAL  LAMED,  MCCURS,  MXT,  MXM,  SDIST,  LAMBDG 

INTEGER  SGNCGT 


I  EN    =-1 

I  EP    =  1. 

I  RSPEED  =  MCD  (  (JBUN-1)  ,15)  +1 

|  KCOURS  *IfII  (JSUN-1)/15.)  +1 


I 

•  CA1L  GGOE  <ISEZD2,2,U) 


CALC0LAI2  ESTIMATE  OF  TABGET  C0URS2  (ST2P  31) 


19 


|  IACM  =TAC+C2 (KCOURS) 


20  |  EIFCO=TACM-IAC 

I 
(CONTINUED  ON  PAGE    2) 
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?AG2  2 


.     *         *    . 
♦If  * 

TAC21.GE.EB2 

*    .  .    * 


.     ♦         *    . 

*  IF  * 

TACH.II.O. 


"^  i     *p 


I    TACa=TACa-PH2 


zf, 


'  -^ 


*   .  * 
I 

E 


.  » 


J  TAC^=PH2fTACa 


I 


CALCULATE  EST  Of  TA&G2T  SPEiD  (ST2P  32) 

_-L. 

|  TA.1   =  TA>S2  (X5PE2D)  

I 


CALCULATE  TOSFECC  CEfLZCTICN  AHGL2  (ST2?  33) 


I  ASP   =EEA3-:ifCO 


.  *    *  . 

.  *     If      *  . 

*     ASP.ll.-EHI     * 


I  ASP=?H2+ASP 
I 


IACS  *IAH*SIH(ABS(1S*) 

IAA   =TACS/TC 


(CONTINUED  CN  PAGE    3) 
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I 

.  *     II      *  . 

*  EAA.GE.1.      * 

*  .  .  * 

*  .        .  *          I     I  26    | 
*  .  * 


j  EA    =  A£SIN(EAA) 
I  £A    =  SIGS  <EA,ASE) 


CALCOLA12  TCSPEDC  JJAIH  FIBING  COURSE  (STEP  3U) 
__[_ 

|  ECOUHS  =EHC*EHI*DA  | 


10  -  *    *  . 

.  »     IE      *  . 
*    aC0aES.Gl.PH2    * 

*  .  .  *  r         i   scouas=acoofls-?H2 

*  .   *  


«   *       *  . 
-   *         if  *  . 

*  i3CCUB5.Ll.0-  * 

*  .        .  «  I     |  SCOUSS=PH2+HCCOHS 

*  .  *  

I  I 

i 
.  *   *  . 

.  *   ii    *  . 

*  MC0UBS.G1.PH2    * 

*  .  .  * 

*  .        .  *  I     |  10    | 

*  .  *  


(CONTINUED  ON  PAGE    H) 
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CALCULATE  TOEEEDC  PRES ENT. FIRING  COURSE  (STEP  35) 
I 


I  CIFALF  =  <U  <1)  *2.*AIFA) -AL?A  | 

i  TCOUHS  =flCCOES*DIFALF  | 


.  *    *  . 

.  *     IF      *  . 

*  1C00BS.GI.SH2    * 

*  .  .  *■ 


*  .        .  *  r     |  TC0URS-TC0UES-PH2 


.  *    *  . 

.  *     II      *  . 

*    ICOORS.12.0.     * 

*  .        .  *  T     |  TC0URS=PH2+TCCURS 


I 


.  *  !  *  . 
.  *     II      *  . 

*    DIPALI.GE.C.     * 

*  .        .  *  T     |  EN=-PN 

*  .  *  , 

J 


CALCULATE  ESTIMATE  0?  TARGET  RANGE  (STEP  36) 


.  *    *  . 

.  *     II      *  . 

*  0(2)  .GE.C.5     * 

*  .  .  * 

*  .        .  *  P     |  PP=-PP 


*  .  * 


I 

I 


|  5HGDI?  =  (1  .-U(2)  )  *RANG2*0.  15 
(CCNTIHDEE  ON  PAGE    5) 
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HKGHOD=BANGE+SIGN(BNGDIF,PP) 


END 


IACMG='IiCa/5AD  | 

TAM1  =TAH*2  | 

IA1   =C£/HAC  I 
EA2   =Ci2*CA1 


PRINT  CUT  OF  FIRING  DATA 


.  *    *  . 

.  *     IF      *■  . 
*     IP2INI.SC..1     * 

*  .        .  *          T     |  25    | 
*  .  *  

J 

I 
I 

***WRITZ  (6,  122) TAC3G,IAai  ,ZNGaOD 

122  FOBHAIMIf'EST    OF    TARGET    DATA    FOR    FIHING* ,/, 

4X, 'COURSE'  ,5X, 'SPEED'  ,61, 'RANGE'  ,/, 
11,3  (F8.  1,35) ) 

♦♦•WHITE  (6,  12a) CA1 

12U  FCJflAT(U,'TCBE    DEFLECTION    ANGLE    IS  «,F6.2) 

_i_ 

I     COOH    =SCCCES/RAD  I 


*»*WRIIZ {6,  125) COOH 

125  FCflaAT(U, 'TORPEDO    SAIN    COURSE' ,8X,F6.  2) 

126  CCNIINOZ 

I 

I     25   ~T 


26  «»*WRITE  (6,3C) 

30                         FCEaAT(/,1I,'NCT    FEASIrLE    TO    FIRE    DOE    TO    LOH    TORPEDO    SPEED') 
j_ 

|     IFLAG=1  | 

I 

25  BE10RN 
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I 
A  TCEEEDO  SIMULATION.    SUBROUTINE  POSIS. 

SOEEOUTISE  ECSIS 
IS  CALCULATING  NEM  POSITIONS  OF  TAfiGET  AND  TORPEDO  IN 

EACH  llfli  STEE 


COHMON    ISEEC2,    TTiaS,     TO,    TA ,    TRATE,     RANGE,     ALFA,     LAKSE,    TADEC, 
3EAE,    1AL,    IAC,    CCC3,     DEVSP,     BNG,    PN,     ?H2,     OCCURS,     TCOU3S, 
3il,    3Xa,    IK,    IDTI2E.    ITIHS,     XT,    it,    xtar,     TTAfi,    TDIST, 
SEISI,    TUENTC,    INTVAi.,    PHI,     3MAX(5,o),    TRANGE,    0I5T,IPSINT 

REAL    LAHEE,     HCCUHS,     2XT,    (US,    SDIST,    LAaBDG 

INTEGER    5UNCDT 


TIME    CCUNT     (STEP    C1) 
_]_ 

|     IK         =IK+1 

I     IIiaE=ITiaE+IETIHE 


CALCULATE    TOTAL    TCSP    3UN       AUD    TABGET    HON     (STEP    C2) 


I     3XT      =C1T*1£IST 

I   tin     =exa>aDisT 


CALCUIAIE    NEW    ECSITION3     (STEP    C3) 
I 


|     XT  =XT  +  SIN  (TCOUR.S)  *TDIST 

IT         =YI*CCSJTCOUE3)  *TDIST 
|     XIAR    =XIAR*SIN  (TAC)  *ttDISr 
{    'iTAR    =  ITAH*CCS  (lACi  *MDI3T 

1 


CALCULATE    SEH    TOSP    COUHSE     (STEP    C4) 


J     IXCOUR    =TCC0aS  +  5IGJJ  (TUaUTO,EN) 

I    IX-CDIF   =A£S  (acouss-Txcoua) 


(CONTINUED  ON  PAGE    2) 
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.  *    *  . 

*     If      * 
TXCDIF.GT.PHI 


|  TXCDIP=PH2-TXCDIF 


ENE 


.  *    *  . 
TXCDIF.Ii. JIFA 


I  15    I 


I 


I  EN    =-£H 

3LPADI  =TXCDI?-AL?A 
t  1XCOOH  =TXCCUR+2.*SIGU(AL?ACI,PN) 


15 


|  1CCU2S  =TXCCUH 


.  *    *  . 

*     If      * 

TCOUB5.GI.PH2 


*  .  * 
I 
F 


|  TCO0aS=TCO0ES-?H2 
I 


25 


.  *    II      * 

»   icooas.ii.o. 

*  .      .  * 
*  .  * 

I 

F 


HZ10RN 


I  TC0U2S=PH2+TC0URS 
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\ 
A  TOBFEDO  3ISULA1ICN.    SUEHOOTINE  DETECT. 

SUEBOUTINE  DETECT 
TO  CHICK  IF  T4BGZI  IS  DETECTED  AND  STOBZ  DETECTION  DATA 


COMMON  ISZED2,  TTIME-,  TO,  TA,  TRATE,  BANGS,  ALFA,  LAMED,  TADEC, 
EEAB,  BAD,  TAC,  CCCK,  D2V3P,  3NG,  EN,  ?H2  ,  MCCDBS,  TCCU3S, 
MIT,  MXM.  IK,  IETIME.,  ITIMS,  XT,  KT.  XTAB.  YTAS,  TDISI. 
HEIST,  TUENTC,  INTVA1,  PHI,  BMAX(5,o),  TBAllGE,  DI3T,IPEINT 

CCM MO N/DA1 A/ El,OFLCE 

BEAL  LAMED,  flCOUBS,  MXT,  MXM,  MDIST,  LAM3DG 

INTEGEB  BUNOUT 

DIMENSION  EE  (3,2) 

DOUBLE  PBECISICN  PC*MA2,B , BBX,BEX1 ,BBX2, X1 , X2 ,X3,XX 1 , 

V1,V2, V,U,FIf ACI,XC1,XQ2,X0  3 ,XX2 , ?OW ZB ,2  2L3 , ALAM 

SETTING  OF  TABGEl  DIMENSION,  A  -  IABGET  LENGTH, 

E  -  TAEGET  WILTE,C   -  TAHGZT  DEPTH.  (STEP  D1) 

|  A     =100.  I 

i  E     =15.  I 

I  C     =4.  | 


CHECK  II  TRANSMISSION  (STEP  D2) 

\ 

I 

.  *    *  . 

.  *     IF      *  . 

*    IK.LT.INTVAL     * 

*  .        .  *  T     \.    20 

*         *  


I  IK    =0 
(CONTINOID  ON  PAGE 
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I 

CALCOIAIE  SANC-E  TC  TARGET  (STEP  D3) 
I 


EIFX  =X1AB-XT 
DIFY  =YIAR-'?T 
LIST  =SCBT  (  <DIFX**2)  ♦(DIPY**2)  ) 

I 


TEST  If  TARGET  IS  BIIHIN  POSSIBLE  DETECTION  RANGE 

( 

I 

.  *    *  .  ' 

.  *     IF      *  . 

♦DI21.GT.  (TAEEC+A/2.) 

*  .  .  * 

*  .        .  *  T     |  20    j 

i 

F 


DETECTION  THRESHOLD  (STEP  24) 
I 


PCSJMAX  =  1./  (3*TAD2C)  **« 
U     =SCAL£  (PHI/2.) 
CONST=  (A**2)  *  (B**2)  *  (C**2) 
POHdAX  =CCNST*POSaAx«U 


CALCULATE  TIUE  TC  TARGET  (STEP  D5) 
I  3IHDL.1  =0121/1500. 


I  liaDLI  =0121/1500. 
»  XTAR1=X145+SIN  (XAQ  *  (TIMDL1*TA) 
YTAR1=Y1AR+CCS  (TJC)  *  (II3DL1»TA) 


CENTE  -  HEARING  CF  CSNTZR  OF  LOBE  (STEP  D6) 


ED    =-CA»CFLOB 
CENT3=TCC0ES+CD 


I 

(CONTINDED  ON  PAGE    3) 
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I 


PAGE 


.  *    *  . 
.  *     IF      *  . 

*  C2NTE.II.0.     * 

*  .  .  *- 

*  -        .  * 
*  .  * 


|  CENTB=PH2+CENTB 


t 

I 
.  *    *  . 

.  *     IP      * 

*    CENTS. G1.EE2 


|  CENTB=CENTS-PH2 


CALC0IAI2  EEAFISGS  TO  TA3GET  (STEP  D7) 


EB    =  EEAFIH  !XTAS1,ZTAE1,XT,YT) 
XTAB2=XIAfl  1+SIN  <TAC)  *A/2. 
YTAH2=YTARH-COS (TAC$ *A/2. 
CISI1=£CSI  I  (S1AE2-XT) **2> (YTAB2 


HB1 


-IT)  **2) 

=EEAEIN  <XTAB2,*TA32,XT,7T) 

__  e  ■•  \t    i  m  ■»  A  *  *  *  /"• 


XTAH3=XTAB  1-SIH (TAC)  *A/2 
'iTAa3=YTAS1-COS  (TAC)  *A/2» 
EIST2=SCHT (  (XTAE3-XT) **2*  (TIAH3 

-*T)**2f 
SB2       "E2ABIH (XTA33 , ZTAH3,2T,IT) 
EEL1    =HE1-CENT3 

I 


*  II  * 

azn  .gi.ehi 


I     R£L1=BEL1-?H2 


(CCNTINOED  ON  PAGE    4) 
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A  TOBPEDO  SIMULATION.    SOEIOOTINS  DETECT. 


PAGE    4 


I 
.  *    *  . 

*     IE     * 

REL1.LI.-EHI 


*  .  * 

I 
f 


T     |  H2L1=PH2+HEL1 


|  5EI0  =E£-CESTB 


.  *    * 
*     If    '  * 
RELO.GI.PHI 


T     |  RELQ=R2L0-PH2 


.     *         *    . 

*  IF  * 
3ELQ.ll.-tHI 

*  •  .     * 


I     HEL0=PH2+S2L0 

I 
I 


I  5EL2  =  HE2-CENTB 


.  *    *  . 

.  *     I? 

*     HxL2.GI.PHI 


(C 


1 


CHTIUDZD  OH  PAGE    5) 


T     |  HEL2=HEL2-PH2 


Hi8 


A  TOEFEDO  SIMULATION.    SQEECUTINZ  DETECT.  PAGE    5 

I 

I 
.  *    *  . 

.  *     If      *  . 
*    Bfl2.lI.-EBI    *   

*  .        .  *  T     |  R2L0=PH2*-R2L2 

*  .,  *  - 


CALCULATE  TBANSHISSION  GAIN  FACTOR  (STEP  D8) 
1         

i   aaxi  =  beli 

I     B3X      =3110 
I    SBX2    =  B2L2 

I   iLia  =iah2D 


I 


COME01E    SEFAEATE    GAIN     JACTCBS     (STEP    D9) 
I 


X1  =  XF  ACI  (33X1  ,A1AK} 

X2         =XrACl  iHBX.lLAa) 
13         =XEAC1  jEBXi,ALAa) 


XXI       =  (X1  +  X2+X3)  /3.  | 


CALCULATE    TARGET    ASPECT    AND    TARGET    SCNAB    CHOS5- 
A32A     (T'RGET    STRENGTH)     DOE    TO    ASPECT     (STEP    D10) 


I     ANTB3=5E*FBI  | 


.  * 

.  *     IF      *  . 

*    ANT3B.GT.FB2     * 

*  .         .  *  r     |  ANTBB=ANTBB-FH2 


*  .  * 

F 

I- 


|     EELA    =ASTEE-IAC 
(CONTINUED    ON    iIgE         6) 
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A  T0EE2E0  SIBULAIICN.    SCEE0UTIN2  DETECT. 


PAGE 


.  *    *  . 

*     IF      * 
BELA.GT.EHI 


I  H2LA=PH2-3ELA 


.   *         IP          *  . 

BELA.ll.-JHI 

*    .    * 

: 

i 

I  BELA=PH2+BZLA 


I  EELE  =  AES  (EEIA) 

V1    =  (A**2)  *  (DCCS(5EL3)  **2) 
V2    =  (E**2[*  (DSIM.  SEL3)  »*2 

i  7     =  (V1+V2) **2 

I 


COaPOTE  SCALING  JACTOB  DUE  TO  ASPECT  (STEP  D1  1 ) 


0     =SCALZ(BEL3) 
FIFACT  =0/V 


CALCULATE  BEIUBNTISE  ?CR  ECHO  (STEP  D12) 
I  TIHDL2  =2.'TIMDL1 


CALCULATE  HEL  EEABIMG  FOB  H2T0H3ING  ECHO  (STEP  D13) 
_J_ 

I  TOBHST  =TBATE*TiaDL2  | 

|  TXC   =TC0DSS+SIG;J  (TUENST,?H)  I 

TXDC  =A£S  (MCCUES-OC)  | 


(CCHTIHUEC  CN  FAGS    7) 
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A    TOEPEDO    SIflOLJTION.  SDEBOOIIUB    DETECT. 

I 


PAG2         7 


.     *         *    . 

.     *  IF  *    . 

*  T2DC.GT.E8I  * 

'     *    .  .    *  T 

*    .     * 


I     TXDC=PH2-TXDC 


18 


19 


*     IE      * 
TIDC.II.AIFA 


I  18    | 


|  ZLF      =TX£C-4LEA 

I  S£LTO=SEL0-2IGM (  (TIJRNST- ALE)  ,?N) 

I        +SIGS  (ALF,EN) 


I  19   I 


i  EELTO=F£ia-SlGH(TU-RNST,PM) 
I 


SELDIZ  =BZII0-RSIO 

E£LT1=BEL1*BE1DIZ 

8ELT2=EEL2+5ZLDIZ 


I 
.  *    *  . 

.  *     II      * 

*  PELIO.Gl.Efll 

*  . 

I 


T     |  aEIT0=PH2-BZII0 


(CCHTINOED  OH  EAGZ 


3) 


151 


A    TOEPEDO    SIBULATION.  SCEflOOTINE    D2T2CT. 


?AG2         8 


.     *         *    . 

*  IF  * 

RELIO.II.-EHI 


•     *         *   , 
*  IF  * 

BELI1.fil.EHI 


I     RELTO  =  PH2-t-3SLT0 


|     H2LT1=PH2-BELT1 


*  IF  * 

E2L11 .Ll.-EHI 


I     2ELT1=P»2>RELT1 

I 
I 


.  *    *  . 
*     IF      * 
HELT2.G1.EEI 


T     |  HELT2=PH2-BELT2 


.  *  '  *  . 
*     IF      * 

HELI2. II. -E8I 


I  RELT2=PH2+aZLI2 


(CCHIINOED  C 


S  EAGE 


9) 
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A  TOEPEDO  SIMULATION. 


SQEBOUTINE  DETECT. 


PAGE 


I 


11S1    FCE  DOPFIZE 


I 


F3QDEF  =50.*2.*(TO*CCS  (CD)-1.) 

/150C. 
ESQDI?  =50. *2. *(TO*CCS(DD) +1.) 

/1500. 
EELC  =  E£LE 
I  FEQSH=50.*2.*{TO*COS  (DD>  +TA*COS  (RELC) ) 

/150Q. 


.  *    *  . 
.  *     I?      *  . 

*  (FEQSH.LI.FBCDIF)  .AND.  (FEQSH-GT. FECDEF) 

a  m  m       *__ ___ — 

*  .        .  *  T     |  20    | 


CALCULATE  3.ANGE  AND  3EA3ING  TO  CLOSEST  PART  OF  TABGET  (STEP  D1U) 
l_ 

=DIST  I 

=HEIT0+DD  | 

=DISI1 

=BcII1+DD 

=DIST2 

=BZLT2+DG 


I    ca  < 

1,1) 

C3 

1,2 

ZB  | 

2,1 

tB 

2,2 

I     IE| 

3,1) 

11 

3'2l 

zn 

$    =C 

1     3D 

=  1 

I 
.  *   *  . 

.  *    If    * 

*   C3XH.GI.EB  (2,.1) 


10 


I 

I 

.  *    *  . 

.  *     IF      *  . 

*  CaiN.GI.EB  (3,  1)   * 

*  .  .  *- 

*  .        .  * 

*  .  * 

I 

F 
I 
I 

(CCNTIN052  OH  PAGE   10) 
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A  T0BP2D0  SIHaLlIICH.    5CESO0TINE  DET2CT. 


?AG2   10 


I  12   | 


10 


|  DHIS  =  CE  (2,1) 
J  3D    =2 


I  9    | 


11 


12 


CaiS-  =  CE  (2,1) 
fiD    =3 


I 


|  EBEI  =CE  (.1Ey2) 


.  *    *  . 

*     II      ' 
BHEL.G1.5BI 


.  *    *  . 
*     If      * 
BSZI.I1.-EHI 


I  BBEL=PH2-3BE1 


I  BBEI=PH2+BHEL 


CALCOIATE  3ZCZIVING  GAIN  iACTOH  (ST2P  D15) 


I  EBX1  =HEIT1 
|  EBX  =SELTC 
I  BBX2  =SEIT2 


X01  =XEACI  (BEX1  -ALAS) 

I  X02  =XEAC1  (B3X.ALAH) 

|  102  =11  ACT  (R3X2.ALAM) 

I  XX2  *(IQ1*1Q2+XQ3)  /3. 


(CCH1IHDID  OS  EiGE   11) 
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A    TOBEEDO    SlflOLATICN.  SCEBOUTINE    DETECT.  PAGE       11 

I 

CALCULATE    FRACTION    OP    POiJEB    IN    TO    RECIEVER     (STEP    D16) 
I  

I    E0HEB=CCNST*IS1*XX2*FIFACT/ (DI5T  I 


TEST    FCH    DETECTION    THRESHOLD     (STEP    D17) 


.    *         *    . 

.     *  IF  *    . 

*  EOHER.  LI.ECS2AX       * 

*  .  .     * 

*    .  .    *  T  |     15         | 


COMEOIE  HEARING  FATE  (STEP  D13) 
I  ASP   =8E+PHI-TAC 


.  *    *  . 

.  *     IF      *  . 

*     ASE.G1.PBI      * 

'  *  .        .  *  T     |  ASP=ASP-?H2 


*  .  * 


.  *    *  . 

.  *     IF      *  . 

*  ASP.L1.-EHI     * 

*  .  .  * 

*  .        .  *          r     I  ASP-PH2+AS? 
*  .  *  


| 


|     AP         =ASP*fiELIQ 
(CCNTINOED    ON    PAGE       12) 
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A  TOBPEEO  SIEULATION.    SCZBOUTINE  DETECT.  PAGE   12 

TAC£=TA*SIN  (AES  (ASP) ) 

J  IOCS  =TC*SIN  (ABS  (EELTO) )  | 

EaATE=  (TACS+SIGN  (TOCS,AP) )/DIST        f 

I 
CHECK  SEARING  HATE  AGAINST  TOHNSATE  (ST2P  E19) 

I 

.  *    *  . 
.  *     II      *  . 

*   EHATE.GE.I3A12    * 

*  .        .  *  I     |  15    | 

*  .  *  


CHECK  TOBPEDC  SPEZC  ADVA-NTAGE  (STEP  320) 
_J_  

|  TAIS  =T**CCS  (ABS  (ASP)) 
|  10LS  =TC*CCS  (ABS  (RBI1D)  ) 


J 

.  *     IE      *  . 

*AES  (ASP)  .Gl. PHI/2.  * 

*  .  .  * 

*  .        .  *  T     |  TALS=-TALS 


I 

I 

.  *    *  . 
.  *     IE      *  . 
*AES (2ZLTC)  .G1.2HI/2. 
*  .  .  * 

*  .        .  *  T     |  TOLS=-TOLS 


(CONTINUED  ON  EAGE   13) 
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A  TOHPEDO  SIMULATION.    SCEHOUTINil  DETECT. 


-  *    *  . 

.  *     IF      *  . 

*  (1ALS+TCIS)  .1E.0.  * 

*  .  .  *- 

*  .        .  * 


I  15    | 


PAGE   13 


HI£  =BE-accoes 


I 
.  *   *  . 

*     IF      * 
315. GI. PHI 


I  HLB=aL5-PH2 


.  *    *  . 

*  IF      * 
HLB.II.-EHI 

*  .        .  ' 

*  .  * 
I 
F 


RLB=PH2*ELE 


SIOBE  CSTECTICN  £ATA  (STEP  D21) 


JCONT=JCCNT*1 

jaAX    =  2  310  jJHAX,  JCQNT) 

ICONT  =  JCCiNX 


(CCJJTINUFE    GN    PAGE       1ft) 
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4  TOEPEDO  SlfiOLATION.    SCEHOUTINE  DETECT. 

I 


PAGE   14 


.  *    *  . 

.  *     IF      * 

*  ICOH1.GI.5 

*  . 

*  .        .  * 


I  ICONT=5 


25 
30 


STOEE  DATA  IN  ACCCEEANCE  WITH  SUAEER  S0CCSSSI7E 
DETECTIONS  (S1EP  D22) 

GO  TO (30, 31, 22, 33, 3U)  ,ICONT 

I 

.  *    *  . 
.  *     If      *  . 

*   ESAX  f1  ,1)  .KE.O.   * 

*  .        .  *  T     |  20    | 


I  FHAX(1,1)  =CIST 

|  5aAX(1/2)  =  (RELTC+DDJ/RAD 


.  *    * 
.  *     IF      *  . 

*  RflAX  (1 ,2) .GI.  180.  * 


*  . 


I  5HAX(1,2)  =C3IH 

BfliX  1,4)  =8EEL/5AD 

I  EMAX(1,5  =RELA/RAD 

I  SHAX(1#6)  =EIE/RAD 


I 

20 


I  HMAX  (1  ,2)  =  2MAX  (1,2)  -360. 


(CCHTINOZD  GN  E 


Age 


15) 
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A  TOBPEEO  SIKOLATICN.    SCESOUTINB  DETECT.  PAGE   15 


.  *    *  . 
.  *     IF      *  . 

*   rflAX  (2,  1)  .  IZ.O.   * 

*  .        .  *  T     i  20 

*  .  *  


EH1X(2,1)  =DIST 

BaAX(2,2)  =  (HELT0-K)D)/RAD 


.  *    *  . 

.  *     IP      *  . 
*  3MAX  (2, 2) .G1. 130.  * 


J  SMAX  (2,2) =BHAX(2,2)  -36C. 


BaAX(2,3)  =CKIH 

E3AX(2,4  =E5EL/HAD 

SaAX(2,S)  =HELA/aAD 

SaAX(2,6)  =5L£/3AD 


I  2C   | 


32  .  *    *  . 

.  *     IF      * 
*   SHAX  (3,  1)  .HI.O.   * 


4 — 

*  T     I  20    | 


J 

I 

I 

_i_ 

|  5MAX(3,1)  =DI3T 

RKAi(3,2)  =  (3ELT0*0a)/HAD 

(CONTINUED  ON  PAGE   16) 
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A  T0EP2D0  SICOLAXICa. 


SCEEOUTINE  GETECT-. 


PAGE   16 


33 


.  *    *  . 

.  *     If      *  . 

*  E2AX  {3,2)  .G1.180. 


I  H3AX  (3,2)  =HflAX  (3,2) -360. 

_ 


I 


EBAX  (3,3 
EaAX  (3  ,4 
EMAX  (3,5 
EaAX  (3,6 


=caiN 

=  EEEL/SAi> 
=BELA/HAD 
=  ELE/'aAD 


I  20   | 


.  *     IF       * 
*   S3AX  (4,1)  .  KE.O. 


I  20 


BSAX  (4,1) 
EHAX  (4,2) 


=  CIST 

=  (RELTO+QDJ/EAD 


.  *     II 
*  SfiAX  (4,2)  . 


*  . 

1 

r 


G1.180.  * 

.  *- 


EMAX  (4,2)  =33 AX  (4,2)  -36  0. 


EBAX(4,3) 

=  caiN 

i 

E3AX  (4,4 

=  EBEL/?AJ3 

j 

EaAX  (4,5 

=3ELA/3AD 

=6L£/BAD 

I 

fiaAX  (4,6) 

1 

(CONTINUED  ON  PAGE   17) 
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A    TOBESDO    SIEOIAIICtf.  SCEECOTIIiE    DETECT.  PAGE       17 


J     20       i 


34  .     *         *    . 

.    *  IF  *    . 

*       E8AX  (5,1)  .  SE.O.       * 

*    .  .    *  T  |     20         | 


B3AX(5,1)     =DIST 

BaAX(5,2)     =  (flELTOfDD)/RAD 


I 
.     *         *    . 

.     *  I?  *    . 

*  EiJAX  (5,2). GI. 180.     * 

*  .  .  •  * 

*    .  .     *  ?•  |     33AX  (5,2)  =SHAX  (5,2) -360. 


*    .     * 


i 


i 


i  EdAX 

(5,31 

=oaiN 

EMAX 

b'u 

=EEEL/HAD 

j  EMAX 

5,5 

=BELA/BAD 

1  ESAX 

b,6 

=SLE/RAD 

I    20       i 


IF    NO    DETECTION,    SET    DETECTION    5TAI0S 


15  |     ICQNT=Q 

I     JCONT=0 


20  B11UES 

END 
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A  TOHPEDO  SI8ULAIION.    FUNCTION  3ZAHIN.  PAG2    1 

I 
A  TOHPEDO  SiaOLAIICS.    FUNCTION  2EARIN. 

FUNCTION  cZAHIN  (A,3,C,  D) 
TO  CAICOLATE  SEABING  FHCfl  I03PSD0  TO  TAfiGET 
I 


EIFX  =A-C 

CIFY  =3-D 

EH2   =  2. *i.  141592654 

BAD   =PE2/260. 


1 

.  *    *  . 

.  *     I?      *  . 

*  DIFY.NE.O.      * 

*  .  .  * 

*  .        -  *  T     |  16    | 


|  HE    =90.*5AD 


I 
.  *    *  . 

.  *     IE      *  . 

*  DIFX.LI.C.      * 

*  .  .  * 

*  .        .  *  I     |  3B=fifl+  (180.  »3AD) 


*  .  * 

i 

I- 


i    n     i 
_JL 

16  |  HS    =AIAN2  (3I?X,DIFY) 


(CONTINUED  ON  PAGE    2) 
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A    TOEEEDO    SIMULATION.  rQNCTIOH    32ARIN.  PAGE  2 

I 

I 
.    *         *    . 

.     *  IF  *    . 

*  B8.I1.0.  * 

*  .  .     * 

*    .  .    *  II     HB=R3+PH2 

V ~'l 

-I- 

17  |     EIAHIH    =HE  ! 

I 

SZ1DBN 

END 
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A  TORPEDO  SICOLATICN.    FCKCTION  XFACT.  PAGE 

I 
A  TORPEDO  SIMULATION.    FUNCTION  XPACI. 

FONCTION  If  ACT  (X,  I) 
CALCULATE  3EE0CTICN-F ACTOR  IN  TRANSDUCES  GAIN  DUE 

TO  RELATIVE  EEAEING  OFF  CENTER-HEADING  CF  TORPEDO 

DOCELE  PRECISION  X  ,  XFACT,  XT ,  X 
I 


ENE 


J 

EHI 

=  3. 

1<41 

SS2654 



* 

.    * 

IE 
X.5C. 

* 
.0. 

*    . 

*  ■ 

» 

*    . 

1 

r 

* 

.    * 

T 

I  10 


I  XI    =X/'f 

|  XEACT=EA5S  (  (DCOS  (X*0.5)  *DSIN  (XT 
•EHI) )/  (XI'PBI)) 


EE1DBN 

I 


10  |  XEACT=1. 

I 
RETURN 


l6Li 


A  TOBEEEO  SIMULATION.    FUNCTION  SCALE.  ?AG2 

I 
A  TCfiPEEO  SIBOLAIICN.    FUNCTION  SCALE. 

FOKCTICN  SCAIE  (I) 
CALC0LAT2  SCALING  FACTOR  IN  THE  PROCESS  0? 

COHE01ING  TAEGE1  STSENGTH 

DC0EL2  PBECISICN  SCALE, REL3, Z, Y 

I  FHI   =3.151SS265U  | 

I 

I 

I 

.    *         *    . 

.     *  IE  *    . 

*  Y.5T.tHI/2.  * 

*    .  .     *  T  |     Y=PHI-Y 


ENC 


I    Z  =  C25  1635*(Y**2)-0.  18555*Y  I 

I     Z  =Z  +  0.C365*DSia (3.* (Y  +  0. 17U53)  )  | 

+  0.C  15*  (1**2.)  *DSI8  (9.«Y/2.) 

|    SCALE=1./Z  | 


REIORN 
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APPENDIX  C 


DETAILED  RUN  PRINTOUT 


T\PTir\L    3!TL*ATIPr!   "AW    n.TI'T 

pi?,i,"G  ".v:r.E      attack  vigle      t.v.get  Prn,",cc    tatet  s^ee" 

3033.0  -CO.O  270.0  13.0 

■j-nnp^p^    PAfWETEPC 

T<r(j;>   n""T,  °V!GE      Tn.V!5. 1 f!T. 7AI       T^r1"    3Prrr         0''cEn    VC!  Z 
750.03  1.00  iiC.0  3*3. * 

|_3~E   ,,I0T'!  TLT.?!    PATE 

20.0  CO 

T'^PPETICL   'MPTM   nr    TACTICAL    St.'EE^-LV'C         lltfl.I 

txec.eticai.  r^r^nc  ?atio      o.osoo 


5>vi,,n   >«\j;;    l^^C    OPT-OGT    r?nf    CE'TE0    1CAP.  I"G       0.0       7!"P3    n~rlZC7!'" 


•n. 


".I'M  "i;r,r::?.         :       i 

137  Or  TAP.CET  "ATA  rnP  TI!T, 

cmt.se  3peep  v/sce 

25G.3  12.1;  33.*;5.1 

mr.r  nr^L^rTior:  a::cle   13       -12.07 

70P.raEPO   f'M"    C'JRSE  17.33 

".IT!    3TPPPC"    (\?TEr.         13G    SSCVTS 


?X'l    TATA    \S    c'lLLO\!G    AT    Z"^  OF    PX-" 
T"TM   T"PP    11  J!!           2730.0 
PIOT    Tn    TVCET             313.0 

70".?     A-CO^PP                     1^.27'.  .5  Tinn     '/.("nonp 

7^^C;T  r.-coop.p   13771. '5  tapget  y-c"r,',n 

TOP.P   ".Al!!    C"t.'PGE           17.027  7T,f    C",llnGr 
■tO    PETECTIO:!   VAPE    PURIfIG    7!"  13    It.'*' 


l-'.OOO.' 

1  :n»n 


;37.G2: 


EOT    nr    TAnCET    PATA    rOP    "!01"C 
CMJP.CE  3Prrri  n\"CE 

"*  c~     r  *  f.     i  *»  ^n  h     ;, 

TORP    PCrLECTII"    VTl"    10         -l-'i .  C7 
•j-rjnncnn   f4AIJ!    C"(.'?GE  15. 1 3 

P.I7'    OT^P^EH    \?TEr         lfc7    3:C,,f,n 


-ifM    p,'4T«    ,\2    rntLT/G    AT    E"P    nr    V- 
TPTAL    T^PP    P.IT  231,0.0 

niGT    "I    7AP.r.E7  3  01.1 

~o«r  ;;-c,0Pr  i-'r230.7 


T  \r*r*r* 


.fTOP" 


13C77.0 


•ynpp    M\|M    C^I'PGP 


•rtr""»    v  —  pT***ri 
rAnCET    Y-p,,',r>n 
70*"'    CH",SP 


'  1 ft  H    T 

ll.*503 


i".\X!'!l'f'   P2TECT1P''    OVT.EO  VP    !JEAP!"CS 
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